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ABSTRACT 


This thesis Gzaminos the design considerations of RJEi 
Imaging Systems . 

Birst, a system idontif ication probi.em for imaging is 
devoT.opcd in.' terms of the Bloch equations and the rotating 
frame arguments. This shows that 

a) The system identification functions( i.e. the spin 
density distribution and relaxation times T^(r^ and 

T 2 (r^)) are linearly related to the RMR response. 

h) The rf excitation affects the system nonlinearly, 
and appears as an aperture function in the imaging equation. 

The syatem design is shown as a hierarchy, from the 
application level, to the system specifications. The case 
study for a clinical imaging system and a system for study 
of multiphase hio—samples indicates I’esolution of 4*5 mm 
and 1 mm respectively for a frequency resolution of 2 EHz. 

For the latter the gradient and the rf power dissipation in 
the coil arc shown to be 0.5 (r/cm (for a 200 Hz resolution) 
and 250 W pulsed rms and 10¥ G¥ power, respectively. These 
figures are achievablo by using conventional wide-line NMR 
spectrometers. The imaging modality rccommondod is the 
projection reconstruction technique, owing to the ©ase of 
signal acquisition. 
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Eig coupling circuit de-sj^n for tho rf transceiver, 
sliows stringent rGiuiroments in terns of largo power, snail 
post-pulse deringing time (r^} o.5p Sec), impedance Hatching 
with the power amplifier, and receiver SFR. fho use of RIR 
diodes in series if switches, is undesirable due to large 
shot noise ( ^ 100 pA/VHz) introduced. Ihis is one of tho 
reasons for preferring crossed coils in imaging systems, and 
where space is a premium, other configurations need to bo 


oxaninod. 
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CHAPiPBR 1 
HTTEODUCTIOR 

1.1 MR IHAaHTG 

Ruclear Magnetic Resonance laas gained importance as 
an imaging tool, in diemistry eyad Medicine after Lautertur^e 
work in 1973 [^]* H© indicated tliat in the presence of steady 
magnetic field superimposed by a fi^ld gradient, and an 
orthogonal radio frequency magnetic fiold, it is possible to 
extract the information about the spatial, nuclear spin den- 
sity distribution of tho sample under study, fhe method is 
termed Zeugntatography to denote extraction of density infor- 
mation from field interactions (Zeugma means ’that which ;joins 
together’ in G-reek) . Tho clinical applicability of RME was 
shown by Damadian in 1971 [2]. It is seen that certain features 
of HMR make it attractive to clinical imaging, as compared 
to other modalities such as Radiography, Compute d~Axial- 
Tomography (X-ray OT), Uitrasourai, and Positron Emission 
Tomography (PET). A. Moss [3] compares these modalities in 
the light of medical applications. 

Briefly RMR Imaging, is an analytical technique, quite 
similar to RMR spectroscopy. It is able to extract spin 
density infoimation as a function of spatial location, for 
different nude ides such as %(p 2 *oton)^^ isotopes 7£j_ 

13q, 3lp and 27 j* In the clinical context it is usTaally 

tho proton, duo to its inherently large SRR in HMR experiments. 
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In addition there are other ohservable parameters, descrihdng 
the state of the system being imaged. They arc , the spin 
lattice relaxation time, T 2 spin -spin relaxation time and 

V the movement in the sample, all functions of spatial location. 
Through appropriate emphasis on these parameters, different 
aspects of the system such as contrast may be studied (MacovSki 
[4]) . Unlike other imaging methods, spatial resolution in 13MR 
is not primn,rily dependent on the wavelength of the excitation. 
The magnitude of the mcagnetic field gradient, roiorrod hither- 
to, detormines the resolution. 

Tho prop and cons in MR imaging viz-a-viz, other imaging 
tochniiuos is listed helow. Tho basic requiremonts for any 
imaging method in tho clinical context arc given by, Tcr- 
Toggosian [5]« The priiiciplo of HMR laaging is found in 
Oho et ^.,[6] and Kaufman and Crooks [7], while that lor X'-ray GT, 
the technique mainly compared with, is given in Scudder [8] • 

ADVMTAGES ; 

1 • HMR is non invasive and non hazardous (no detrimental effects 
of rf fields in the range of 1-15 MHz, the frequency of 
operation, have been observed) . 

2, The response is chemical specific, enabling tho imping of 
distriihut ion of various chemibal states of the nuclPides. 

5 * It aff ©eta four parameters of moasuromonts, allowing for 
better selectivity of the characteristics to bo measured. 

* vide-^upra 
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4* The dependence of dopt)ler allows selective imaging of 
sample movements, as in flow measurements. This avoids 
motion artifacts found in X-ray GT. 

5» Tlie imaging moda-lities are numerous, with-in the main 

framework of KMR. This introduces flexibility in different 
applications • 

6. The number of moving parts in the MR hardware are fewer 
compared to the high precision linear and rotary motions 
required in Z—ray CT and PET, 

EISADVMTACtBS ; 

1 • The sensitivity of NMR is inherently lower compared with 
other imaging schemes. 

2, The stringent requirements in maintaining high magnetic 
fields (I^5~4 KG) accure^te to 1 in 10^ spatially and 
temporally large magnets for medical applications are 
made with super conductors, which deteimine the economics 
of the imaging system. 

3* In many applications iW'IR is huliy and cumbersome (due to 
the magnet and power supplies) compared to techniques 
like Ultrasound Imaging, which is portable. 

4* Patients with aids such as the Cardiac pacemaker cannot be 
imaged due to the large magnetic fields used. 
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It is predictsd (Eaufmann and Crooks [9]), that while 
Ultrasound and Digital Radiography would continue to occupy 
an important place in diagnostic medicine, IMR would vie 
with X-ray Cl. Though X-rey CT, has been developed and per- 
dBected over a considerable length of time, the commercial 
viability of HMR Imaging with Protons, is imminent in the 
foreseable future. In areas such as in-v ivo metabol in studios 
with 31p and RMR potentially holds the fort|^-. This ai^a 

of research however involves engineering problems (such as 
SHR of the experiment), which needs to be solved. 

1 .2 Purpose of Work 

The ob;ject of this thesis is to present a review of 
different MR Imaging modalities proposed in literature with 
a view of evolving a system design. 

Attention is focussed on mathematical modelling of MR 
imaging. It is compared with the usual linear time invariant 
models encountered in electronic systems. The problem is seen 
in the light of system identification approach in the presence 
of noise-. The relevance of this approach to areas such as 
3 Ip and 13^ imaging is indicated. 

An actual system design for a small imaging instrument 
operating at 3.5 KU (15 MHz) has been carried out. The design 
cons ide rat ioas involve translation of the application require- 
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ments (in this case the study of multiphase b io “specimen) , to 
system parameters. Design of the rf transceiver subsystem 
for this instrument is discussed in deta-il. 

1 .5 ORGMIZATIOK OP THS THESIS 

OEiis thesis considers two aspects of MB. imaging, namely 
the ma.themat ical model, and the system, design. 

Chapter 2, gives the principles of (MR and imaging). 

The va,rious irgpediente in an MB experiment are considered 
and a brief review of different imaging modalities used, is 
given. 

Chapter 3, develops Bloch’s phenomenologica-1 equations 
of MB as the mathematical model. The rotating frame argument 
is employed to derive the basic imaging equa^tion and obser- 
vations regarding the non-linearity of the model, is made. 

The mathematical foimulation for basic excita,t ions, and for 
different modalities, are listed. Finally the processing of 
the MB response, is followed to state the system identifica- 
tion problem. 

Chapter 4, considers, the design considerations for 
the various blocks in the Imaging system. G-eneral design 
requirements are given, followed by the parameter specifica- 
tions for the detection system, consisting of the rf transceiver 
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and the sa^nal aq.aisition system. A case study is presented 
to q.uantify different system req.uirements . The feasibility of 
a small imaging system for use in study of multiphase hio- 
samples is cons idered^ 

Chapter 5, discusses the rf transceiver system, -with 
emphasis on the coupling circuit design and Chapter 6, high- 
lights the problems considered, alongwith some suggestions 
for future work. 

Appendix A presents some details of derivations encoun- 
tered in Chapter 3 of the thesis* 



CJHAPTBR 2 


EB7IS17 OB NMR IMAGING 


2.1 HISTORICAL BACKGROUND 


The Nuclear Magnetic Resonance phenomenon in hulk 
matter was pointed out as ea^rly as la ^936 by, Gorter- However the 
first sucessful KMR experiment was carried out independently 
hy Purcell, Torry and Pound and hy Bloch, Hansen and Packard 
in 1946« Suhseiuently the importance of NMR as an analytical 
tool was realized, and concerted work was underway hy 1950’s 
to develop viable instrumentation. A comprehensive description 
of early IMR spectrometers is found in Andrew [^0]* Chrono- 


logical development in NMR from 1936 till present is given in 


Figure 2.t 



■ II- III 

60 66 70 71 J3 80 82 


Pulsed HMR 
Torrey, Hahn 
1949,1950 

First Successful 
Experiment 
1 946 Purcell et .al • 
Bloch et.al. 

First Experiment 1942 

Gorter 's 
cettemipt 
of'NMR 1936 


Apnl icat ions 
of NMR to 
Medic ine 
Damadian 
1971, 1972 


MR Use of ME 

Imag ing Im^ ing 

Lauterhur in 
1973 clinical 

diagnostic 

1 982 


I : C¥ HMR, Quantum and Classical Theories of HMR 
II: Pulsed MR and Instrumentation 
III: MR Imaging and Applications to C3.inical Diagnostics 


Figure 2,1 A Chronology of HMR 
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The first phase of development was on the phenomenon 
itself, wherein the techniques used involved C¥ excitations* 

The second phase saw the developnent of pulsed techniques and 
Instrumentation. Many of the standard pulsed spectrometers 
were first manufactured during this period. The third phase 
is the application of this technique to imaging and its conse- 
quent use in medicine* 

2.2 THEORY OF MR 

The MR phenomenon arises from the splitting in the 
energy levels of an atom due to the interaction of the nuclear 
spin with an external ma,gnetic field. Heurist ically this 
interaction may he attributed to the magnetic monient arising in 
the nuclei having an intrinsic angulan momentum. This pheno- 
menon essentially occurs in the quantum domain, and must be 
rigorously treated as such. However it is useful to adopt 
the classical arguments to model the phenomenon, especially 
from the point of view of instrumentat ion. For the S8.ke of 
completeness, the quantum Biechenical and the classical pictures 
are presented below [”^0]. 

1 * QUANTUM MSCEANICAI’ INTERPRITATION 

This picture evolves from the concept of the spin angular 
momentum and the spin quantum number I* The magnitude of the 
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angular momentum vector of tlie spinning nucleus is given as 
1 • 

[1(1+1)] ^ where h = h/27i, h being the Tlan]n*s constant. 

!Ehcugh this vector may assume infinite number of directions 
in space, the measurable component (in the direction of an 
externally applied magnetic field) is quantized, as m h where 
m assumes (21+1) values, from I to -I. fhe corx''es ponding 
components of the magnetic moment p(=r e •h/2MpC) of the nucleus, 
are given by mp/I. These components, along with a magnetic 
field , give (21+1) energy levels of magnitude ~ mpE^/I, 
about the main level* This is the energy level splitting, 
referred hitherto. 

The separation between sucessive energy levels 

■ ^o 

^ is, aB = — — 

A radiation of frequency with energy given by 



would lead to an energy level transition mari^d by a net absorp- : 

tion of energy from the radiation. In the case of protons, which 

is of interest in this thesis 1 = 1/2 and hV = 2uH . The 

^00 

splitting of the energy levels and the emission -absorption 
feature are illustrated in Figure 2*2a. ; 
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The discussion thus fa-r has been addressed to the 
interaction of a single nucleus, with a steady external 
magnetic field . On a bulk scale the phenomenon is stati- 
stical and follows the Boltzmann distribution for the popu- 
lation among energj^’ levels. 

Using the concept of the spin temperature^ and the lattice 
temperature (’lattice' refers to the environment, in which 
the nuclear spins are em.bedded) , Time constants, • Spin- 
lattice relaxation time and T 2 » Spin -Spin relaxation time 
and the saturation factor Z, are derived. Bor details the 
reader is referred to Andrew [10]. The parameters T^ and T 2 
would be touched upon again in the next section while discussing 
the classical theory and a . mathsmatical model of MR. 


2, CLASSICAL BTTBRHiSTATIOR 

■■ Hotw ithstand ing the rigour of the quantum mechanical 
treatment, the classical picture gives useful insight into 
the phenomenon. As such this approach is employed in the 
rest of the thesis. A comprehensive trea,tment is given in 
Andrew [10]. 
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In brief, the phenomenon may be Tisualised as follows 

(Figures 2,2 b, c ). 

“ In the absence of a magnetic field the elemental magnetic 
moments (associated with the spin of the nuclei) are ran- 
domly oriented in space, the randomness being governed by the 
temperature of the system, 

- On the application of a magnetic field along the Z’ axis, 
the spins tend to align along the direction of the field, 
the tendency being opposed by random fluctuations imposed 

by theimal agitation, 

- The nuclear magnetic moment experiences a torque proper- 

tional to p 2: and starts process ing about ( or Z 

axis) at the angular rate =^H^, very similar to a 
rotating top or a gyroscope. The constant y, termed the 
gyromagnetic ratio varies with the kind of nucleus partici- 
patizig in the process and is numerically equal to p/I hj'the 
quantum mechanical counterpart. For the proton = 

42,576 MHz/Tesla- or 4*2576 KHz/G-auss. 

- Introduction of a circularly polarized radio -frequency 
magnetic field (at frequency w = yH^) , in a plan# 
orthogonal to the Z-axis, tends to rotate the magnetic 
moment p^ about its instantaneous axis. The rotation is 
termed Nutation in G-yroscopy pai'lance . This is better 
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visualized and handled mathematically hy adopting a 

rotating frame of reference. This frame, teimed X*Y*Z’ 

rotates about the Z axis a,t frequency in the same sonae 
s 

as the precesion of the nuclear magnet, and is implicitly 
h 

assumed in the rest of the thesis. 

Mathematically the precession frequency also called 
the larmour frequency is expressed as 


^ laboratory frame of reference XTZ* 
In the presence of an rf field of frequency Wq, 

res ^ res 

= ~Y (Hq Bf(t);in XTZ frame. 

where denotes the resultant magnetic field. 

res ^ 

In the r'otating frame 


w. 


res p 


= Y 


r> 




max 


as the effect of is subsumed in the rotating frame 


(w^ ) is the Nutation frequency, and is an important variable 


res' r 
in BMR. 


If the rf “field is applied unifoiml^?" for a duration 1^, 


the angle 0 by which nutates is given by 
</! = iy Ip- 
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0 Is temed the flip angles It will "be seen that 
the ohservation plane in IMR is the ZI/Z’Y’ plane, and conse- 
quently to obtain maximum response, 0 must he -n:/2. 

Another flip angle of importance in pulsed KhR studies 
(where rf-excitation is pulsed) is ^0 = These pulses are 
also called the 90° and 180° pulses for obvious reasons. 

Though the picture presented is for a single nucleus, 
the results may be extended to bulk mattei* by studying the 
statistical behaviour. This is incorpora.ted in Bloch’s 
Phenomenological equation discussed below. 
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which is related to the static susceptibility 

Here H is the numher of spins per unit volume, 

Y the gyromagnetic ratio, 

H^is iPhe steady magnetic field , 

ItheS.pin Quantum number (=^ for protons), 

is i)he sample temperature 
s 

The numerical value of is 3*25 z 10 Am /Tesla, for 
protons in water at 3lO°Z*. 

This is thus a 1TJ14PBD MODEL* 

Addition of factors doscribing the spin— lattice and the 
spin-spin relaxation effects lead to a complete set of pheno- 
menological equations called the Bloch equations. The assump- 
tions ttiade in involving the factors and T 2 are that, ( i) 
the effects are exponential and may he obtained by a linear 
superposition, and ( ii) T^ affects the Z-^omponent of the 
dynamical equations while T 2 affects the X and Y components. 
The rationale for these assumptions is based on observations f 
and is that (i)TIio exponential behaviour is consistent with 
a simple quantum statistical description. ( ii) T^ is related 
to restoration of equilibrium, or return to 2-axis in the 
presence of H^, and T 2 i? rela.ted to the phase information 
loss, in the precession of an ensemble of nuclei, under 
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the influence of local field spi*ead* This does not amount 
to restoration of statistical equilibrium of the overall 
spin “lattice system. 

Owing to this^ and T2 are temed "^he longitudinal and the 
transversal relaxation times^ respectively . The assumptions 
in the Bloch equationshave "been examined by ¥angsness and 
Bloch (1953)* The details are given in Andrew [10]« 


The equation in the labora,tory frame is, 



Y M^xH> 



x_ 


To 


1 


y 





where "^x ^ resultant 

magnetic field, and B^^ , B^^ are the components of a circularly 

polarized magnetic field. H^(t) describes the trajectoiy of 
the net magnetization. A deta,iled discussion of the Bloch 
equation is considered in Chapter 3« 


2.4 PRINCIPLE OP MR MAOnTG 

The concept of imaging using MR follows logically 
from the background in the previous section. 

The laimour equation 

w^ = “ Y C 
O ‘ o 
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inplies, that superposition, of a spatially varying magnetic 
field H^(r^) over amounts to a spatial variation of the 
lamour frequency, 

= -Y (H> + H>(r>)) 

The variation is linear if H^(r^) varies linearly 

over • 

’) may he written as 

= ~ G^(r^) ♦r^ where G-Cr^) is the magnetic field 

gradient • 

In the simple case of a field gradient in the x direction 

w>(x) = -Y (H^-h^ 2 ) , where the gradient is imposed 
over the Eiain magnetic field. 

Alternatively, the equation implies a mapping between 
the spatial coordinates, x and and the frequency w(2!*) . 

The primary response in MR imaging is the distribution of "tlih 
nuclear spin,^viiRll7 of protons* This amounts to an image 
foimation with the spin density distribution denoting the 
intens ity of the image • 

It may he noted that the resolution is detemiined by 
the magnitude of the magnetic field gradient G^( x) rather 
than the wavelength of excitation. 
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MR jmaging offers four parameters to be imaged, namely 
p( 2 :,y,z), spin density 


and 


(3:,y,2)] 


Relaxation times 


a?2 (x,y,z)j 

V , velocity of flov. 


Ihe relaxation times depend on the environment is -whicli 
the spins are embe'''.d.ed« It has been pointed out that in 
clinical imaging ( for the detection of lesions ) relaxation 
specific imaging is more effective than spin density distri- 
bution [g ]• Figure 2 .5b,illustrates the nature of and Tg? 
and the relation of , with the existence of tumours. This 
technique is also sensitive to velocity, the effect being 
similar to the doppler in the Radar Context. The dependence, 
finds use in imaging blood flow (and vessels)* 


2*5 A BASIC MR BXPBRIMBITT 

There are three basic constituents in an MR experiment. 
They are 

1..* A steady Magnetic Field along the Z-axis defining 
the operating frequency w^* 

2** A circularly polarized radio frequency magnetic field, 
in the transversal or x,y plane of frequency w^, providing 
the excitation. 
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3’ The response of the sample under study* 

Two modes of excitations used are the 1* C¥ and 
2# Pulsed fonas. These define the two "basic inodes of experi- 
mentation* They are discussed briefly below : 

C¥ EMR ; Early experiments m WIR were conducted 

in this fom. A linearly polarized rf magnetic field is 

coupled to the sample kept in a. steady magnetic field. One 

of the circularly polarized components interacts with the 

nuclea.r spins to excite them to a higher energy level. The 

coupling of the rf field is achieved with the help of a coil 

(e,g. solenoid), configured along with a transmitter circuit. 

The net absorption of energy, by the sample is detected by a 

a change in the impedence level in the coil. Owing to a gamut 

of interact io 00 in the sample (T^, T2, chemical shift and so on), 

the absorption is spread around the operating frequency w^* 

Consequently , the steady aiagnetic field (H ) or the rf frequency 

0 

Wq is swept to cover the whole spectrum. 

[Phis method however has the disadvantage of long 
observation times, and very slow field sweeps due to the 
inhorcnt insensitivity of the oxporimont. 

2* PCIiSBD MR ; With the advent of narrow spectral 
line studies in spectroscopy, 0¥ MR was found to be inadequate 
due to poor sensitivity. The pulsed method on the other hand 
provides the advantage of simultaneous observation of the 
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tfliole spectrum, using the ME system’s transient response. 

A repeated pulse ercitation-ohservat ion seq.ueace improves the 
sensitivity, or signal to noise ratio over the CW method for 
the same observation interval ( Bmst and Anderson [11])* The 
method however requires a Fourier transformation on the tran- 
sient data to obtain the spectrum. This processing is made 
easy in the light of the present technology* 

The technique is especially useful in high resolution 
spectroscopy where, the magnetic field (H^) stability require- 
ments are stringent. The spectrum is recorded in a short time 
compared to the field fluctuations* Typical values of the 
static magnetic field (H^) and rf field are 1 - 5 KG- in 
imaging {and 15 - 25 ES in high resolution spectroscopy) and 
1 - 10 G- (peak value) respectively. The ma,in field fluctuations 
(both temporal and spatial) are required to be within 1 ppa 

of H . 
o 

The pulsed method is assumed in the rest of the thesis. 
The transient signal recorded is called the Free Induction 
Decay (KID) . This is an exponential envelope (with time' const- 
ant T 2 ) modulated over a sinusoidal carrier of frequency w^. 
However in Imaging, when the field gradients are strong, it is 
seen that the decay follows a behaviour, where K 

is "determined by the gradient magnitude. 
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In many samples ^ fact whicli is used to ottaln 

repeated response a slngl© system rolaxation poriod T-j , 

OThe response is called a Spin Echo and is useful for expediting 
data aquisition in imaging* Tliis is discussed in the next 
section. 

2i6 SPIN ECHO : 

QJhis is a phenomenon characteristic of a pulsed KDffl. 
experiment. The concept is understood as follows » 

The envelope of the FID signal decays with a time 
constant ^-j)* '^^® classical picture this 

corresponds to a de phasing of the ensemble of precessing nuclei, 
as shown in figure 2^h)t Concurrent with the phase loss, the 
magnetization returns to the Z—axis alignment at a rate 
deteimined by the spin lattice relaxation. Due to relatively 
long T|., it is possible to excite the sample to reverse the 
dephasing into a rephasing and thus obtain a repeated FID or 
an echo. In pulsed experiments this rephasing is effected by a 
180° Spin Echo pulse. 

Two spin echo pulse sequences are illustrated in 
5 * j_gQ 3 ^Q 2»4- , They are temed the Hahn and the CPMO sequences, 
and are used in pulsed EMH spectroscopy. [€ ] 

Spin Echo may be used to find Tg., as every rephasing 
with, a net reduction in -t vdri>='ve'wa-l magnetization. 


occurs 



24 


due to spin-s-pin relaxation. In ME. Imaging, this is 
primarily used for effective data acquisition and for obtaining 
l2(^^) sensitive images. Additionally the FID time constant 
I2 imaging experiments is considerably saialler than I!2» 
to increased dephasing obtained by the external gradient. An 
approximate expression is given as 

.where L is the spatial extent of the sample, over which the 
FID is observed. Consequently another method of obtaining spin 
echo, is of footed by using gradient reversal [12]., This is 
indicated in Figure 2 » 4 . Given a gradient G(r^), two points in 
space and r^ are separated in frequency by 

A w = Y &(r^) '(^2 - r^”) 

and in phase byA? 2 f(t) = Awt (^Y^Cr^ *" r^) t) . 

If after time t = Tg , the gradient is reversed, 

A0(t) =A wt = Y(*-G(r2 - r]')(t -f2^0 ([||) for t 2: 

thus leading to a rephasing at time t t= Zlg (0(2!l}g)=o) . 

It will be seen later that gradient reversal is preferred 
over spin echo pulses in imaging, owing to the power require- 
ments on the transmitter. 
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2,7 IMAGING MODAIITISS 

NMR Imaging is marked "by a variety of ways of obtaining^ 
a map of the spin density distribution, as against a fixed 
modxlity in X-ray computed tomography. Ihe reasons for this 
is in the basic imaging equation itself, as will be seen in the 
next chapter. Heurist ically the argument is based on the fact 
that, the field gradient in conjunction with the rf oxcitatioa, 
is capable of isolating the NMR response for a small region 
in space.' Thus the various modalities indicate the way in 
which the sample responses are grouped together. Broadly 
the imaging is classified into Reconstructive and Non-recons- 
truct ive methods. The foimer technique follows directly 
from the applications in CT [ 8 ] , and involves reconstruction 
of the image in its entirity from integral projections. 

This method deals with tho nso of Radon transforms. Non 
reconstructive imaging, on the others hand is a feature 
characteristic of NMR. Another reason for such a classifi- 
cation is that, while Reconstructive schemes are marked by 
complexity of software processing, the latter highlights 
the hardware complexity. 

A brief discussion of the two methods is given followed 
by f igures, illustrating the various techniques. Tables 2.1 ,2.2, 
list important features of each of these techniques. 
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1. ESCOISTRUCTB/S TWAC-IITG : This method is similar 
to its counterpart in CT. The procedure for extraction of 
the image may he understood from the basic imaging equation, 

s(t)= / P(x) e ^ dx. (2*1) 

X • 

which follows intuitively from section 2«4. A more detailed 
derivation is given in Chapter 3 of this thesis. The equation 
indicates a Fourier transformation relationship between 9(3:) 
and s( t) , where 

p(3c) = /T p(x,y,z) dydz. (2..2). 

yz 

Thus s(t) contains a plane “integral, (equation 2.2)*_ Kie 
j&nage is thus obtained by invoking the back projection algorithm, 
to reconstruct a 3~D or 2-D image from projections. This 
necessitates scanning the object space, which is also the phase 
space obtained by f idLd gradients . T‘igure2.5, indicates scanning 
of the Ey phase plane (^3^.= ^ ~ "^^y in a radial 

fashion [I3]- 

The pulse sequences for a 2-D and 3“D reconstructive 
imaging experiment are shown in Eigure 2. 6, 

2. EOiT-EBOONSTEnCTITE IMAGING s All imaging modalities, 
which do not fall in the first category, come under this 
category. The sub-classif ications are 
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a) Point and Line Scan Imaging [14] , [15], 

b) Direct fourier transform imaging [‘16] (EPT) • 

a) Point and line scan imaging are brought about by selectively 
exciting a particular region of the sample i*e* a line or 

a point, at a time, and recording the response. This 
process is repeated to Eiap the whole sample region, This 
procedure is marked by an involved data aquisition, though 
further processing is simplified. The data acquisition is 
also slow, which affects the net image acquisition time, 
an important factor in clinical applications. Figure2.5» 

. shows the scanning in the phase plane for 2~-D imaging, while 

Table 2.2, lists the characteristics of the method thereof . 
A more detailed discussion de referred to Cho et.al., [6 ]• 

b) Direct Fourier Transfom imaging involves a line scanning 

. of the phase plane, somewhat «.kin to a raster. The received 
data is then subject to a direct 2-D or 3"“^ fourier trans- 
form for Image formation.The data is processed in a 2~D or 
3~D foim^ defined by the scanning and is achieved by 
switching the x_,y and z gradients. The method is originally 
due to Kumar, Welti and Ernst [ 16] (KWE), and is a normal * 
mode of imaging in clinical contexts, where some modified 
Direct Fourier Transfom metho(3Bare employed as indicated. 
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TifflLE 2.1 

KSCOIJSTRIJCTIVB IMGBTG 
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3* EEC071RY SEQUENCES ; SBNSITI7E IMAGIiiTG 

Elis seq.uence essentially imposes a condition on the 
repetition time of the read sequence, follo'wed by a modality. 

The rationale is that, subsequent to a 90*^ pulse, applied 
for maximizing the FID response , the spin lattice relaxation 
mechanism sets in, which tends to return the magnetization 
to the 2 axis* As a result, another 90® pulso applied at time 
T^ = T j j of the first ptO.se would evoke a different response. 
Sucessive repetitions lead to a steady state response, dependent ’ 
on the repetition time Optimum repetition time f is 

found to be near T^ * In its basic form, the sequence is teimed 
Saturation Recoveiy (SR), and is illustrated in figure 2*6 * 

The exact dependence on T^ is derived in Chapter 3 of the thesis. 
The method is employed to obtain T^ dependent images, which 
are known to give good results in lesion detection, as shown 
in Figure 2, 5b* 

A modif icat ion of this is the Inversion Recovery (IR) 
sequence shoim in Figure 2 *6. This includes a 180° pulse, 
following the FID read out, amounting to an increased T^ 
dependence, and is seen to give better contrast* Heuristically, 
the dependence is owing to a phase inversion in alternate read 
Sequences* It may be stated that Saturation and Inversion 
recovery sequences are not modalities per se, but are used 
in conjunction with them to elicit the T^ dependence. 



4* SPIN ECHO IMAGING (SE) ; 

Applications involving a strong dependence on spin-spin 
relaxation time T 2 » employ this method* Control over the 
echo time obtains a response dependent on 12 ( 2 :) • 

Ihe response is weighted by a factor will be seen 

in the next chapter. 

5. DOPPLER IMAGING ; 

■ : In clinical conto'at MR has been used for imaging 

blood vessels by exploiting this peculiar dependence. Figure 2,6 
illustrates the phenomenon • A selective excitation of 

a small region of the sample of extent ^x for a duration 

ir 

extracts a response from, the sample due to those nuclei, whose 
net outward drift velocity is less than ax /t^ * Ihus excess 
motion may be filtered out at the detection level, to avoid 
motion artifacts*. 

Nuclei with velocities greater than ax /t may also... 

It 

be imaged by receiving the signal from regions adjascent to 
the excited sample volume by applying an appropriate field 
gradient. A detailed discussion of such selective imaging 
technqiues is given by Mac(xmki [4 ]• fo sum up, the potential 
of NMR imaging lies in the four parameters p , , ^2 

which highlights various aspects of an image. 
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2,8 SXGIMIONS BT 3MAGIIG : 

To aid tlie understanding of tlas figure 2,6 , a 

discussion of some salient excitation pulses^ follows ; 

!, rf HJIiSES ; The rotating rf magnetic field is 
tailored to stimulate either the whole sample or specific 
regions* Consequently there are : 

a- : 90° and 180° Broad Band Pulses 

The spectra of these pulses are broad and are designed 
■fco, oxei1>$ the whole sample in the presence of strong field 
gradients- The power requirements on these pulses are large 
and are therefore used sparingly in iiaaging experimontc . In 
most imaging techniques 90° pulses, initiate a sequence, and 
may he applied before imposing the gradients, thereby 
obviating the necessity for the broad band stimulus. However, 
bihad band 180° pulses are used for pi*oducing spin echo, in 
the middle of a pulse sequence, Somotimes gradient reversal 
is also employed for generating spin echo, though the method 
is constrained hy switching spool* 

b- .;90° and 180° Narrow Band Pulses ; 

These excitations are used to selectively stimulate 
a region of the sample, and are important in non-^reconstructive 
modalities. They are employed to effect slice selections in 
2D-IFT hmaging and in select fro line and point scanning 
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procedures. Tlie choice of pulse shaping is complicated by 
the nonlinear nature of MR, and the problem is addressed to 
in Chapter 3 of the thesis* .It may be noted that in certain 
instances image resolution is governed by pulse selection. 

2) RIBLR GRADIBlTfS : ' Two types of gradients are 
used in Imaging, they are : 

a : Read Out Gradient : This is a constant gradient 
applied in a particular direction, usually along the z axis 
for a. duration T >> purpose is to effect the 

read out of the response, and is used in line scan and direct 
fourier transfom imaging. 

b ; Gradient Rulse ; It alludes to the case where the 
duration of the gradient is shorter than the dephasing time 
l/yG^, 30 that the transverse magnetization remains in phase, 
while being modified marginally by the pulse* Its use is in 
Direct Courier Transform or Fourier Zeu^atography* 



GHAPTSR 3 


MHEMATICAI MODEL OM THE MMR IMA.GIMG SLSTSM 

3*1 KSQUIE3MEINTS OP TUB JIATHEMiTiaAL MEL : NMa AS A 
.. SYSTEM IDENTIPICATIOM PROBLEM 

The use of EMR as an itaaging technique is already a 
decade old, and the method has seen the e"?olution and develop- 
ment of various modalities within the main framework of NBIR, 

A review of these modalities were presented in Chapter 2. 

However with the requireffients of full scale systems, to cater 
to a variety of applications in clinical and chemical contests, 
it becomes necessary to unify, or at least to identify common 
points underlying the various modalities* Any meaningful' 
comparision of the imaging method, from the system point of 
view, would necessitate the development of a mathematical 
mo.del of the processes governing the technique. In MMR a good 
starting point is the Bloch equations, covered in chapter 2 of 
the thesis* However it wi3.1 be seen that the nonlinear nature 
of the equations has been a barrier in developing simple quanti- 
tative models to describe the response, especially to arbitrary 
excitation -pulses. The recent past has seen work in mathe- 
matical description of ITMR imaging, Hoult [17], Mansfield and 
Jfeudsley [18] and Caprihan [ 1 c] • 

Before starting the discussion of the mathematical 
model, the requirements on HMR' as a system identification 
problem is examined. The four parameters in this method are 
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9 (?^), !Ej(r^), 8-5i<3. V, The mapping of the spatial 

: information to the frequency spectrum is obtained by the appli- 
Cation of a magnetic field gradient in any general direction r^ 
(the field itself being aligned -with H^) . Consequently, 
considering the MR system (vith the sample under study) as a 
biLack box, the inputs are, 1. The f ield gradient O^r^) and 
2* The rf ezcitatiouB^ (t) -while the output is the PID s?('^) • 

The main magnetic field defines the operating frequency 

of the experiment, and is not presented as an incut, since 
some form of demodulation or frequency translation is inherently 
assumed. Pigure5,*t^ describes schematically the systems 
approach* The following observations may be made at the outset, 
before commencing an analysis with the Bloch equations i 

They are 

1* The relationship between the input and output is 
involved, due to the nonlinearity of the Bloch equations. 

2* The system is multi-input and single -out put, a fact 
to be borne o«t by the model to be developed. 

3 . Any linearity in the system is not obvious on a 
cursory examination. 

4 . The system response, namely the PID is observed after 
the cessation of the rf pulse. This constraint is imposed by 
practical considerations, namely that the transmitted power is 

ifNifoR’v'ATlCsi'll 




KF 
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large compared with the NMR response, thus swamping out the 
receiver. 


5*2 BIOCH SQUAT 10 ATS (Details in the Rotating Fran© of 
Reference in the Presence of Field Gradients) ; 


The section 2-5 presents the phenomenological model due 
to Bloch* The rationale and assumptions involved therein are 
indicated. The vector dynamical equation of magnetization 
(assuming a lamped model) is given by 


= V 5S2> - ^ 


A 

i 




A 

V 


'f. 


z 




where i + B. (t) i + B. ft) "i^* is the resultant 

o z ts X ly y' 

magnetic field. 


In this equation, the first term x is due to the 
interaction of the net magnetization with the resultant magnetic 
field| while the remaining terms are included phenomenologically. 
These are a set of three scalar equations, detei-miniag the 
trajectory of the x,y and z components of the magnetization fK 


In the case of field gradients being applied to the 
main field, it is more convenient to treat the Bloch equations 
in the rotating frame of reference, as explained in Chapter 2. 
The parameters of the reference frame are as follows t 
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^ * Tlio rotating fr-amo is labelled X*, T* and Z»» 

2 4 IPhe rotation fr^^^quoacy is the Larmour frequency 
given by = yH^* The direction of tbe corresponding vector 
is opposite in sense to H^« 

" 3 .*fnie residual field due to tbe field gradient, at 
location r^ is denoted by AH(r^) . 

4 * fbe rf excitation 3 ^(t) bas a frequency w^- 

Hence the vector equation in tbe rotating frame is 



V 4 ^ 
fa ^ 


LV: 

fa 


■M ) 


A 


(■=^. 2 ) 


where H^* = ^H(r^) i^.' -f- B^(t) i^’ -h 0 .^^ 

. ■ (obtained by choosing the X’-axis along 

the B^(t) phasor) 

and . ^ ^ ^y* * *^z ' » components 

being functions of, time(t) and position r^ (w.r^t^ the laboratory 
frame), siting ^he three components separately yields 

dH » • ' 
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(3.3) 


Def ining: M as 

H = M > - 3 M > 

jr 

(s.ince orthogonal components in the trans'verse plane of the 
rotating coordinates has a quadrature phase relationship in 
the laboratory - frame) , the Bloch equations simplify as , 

(3*4) ; 


The term corresponding to T^ Eiay be neglected for all t << T^ . 
It may be noted that tho Bloch equations are written for a 
spatial location r^ in the sample , so that ^y** * 

T-j ’ and T^’ 9-^® functions of r^ (the 3.ocation in the la.boratory 
frame) • 


+ (|“ -* j = -j yBi(^) 




dt 


= YB^(t)ImM 


iK' “ 




-1 


dB. » 


M.. 


= YB|(t)l!2» - ^ 


dt 


= - YB5.(t)H^’ - 


^ q) 

T. 
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The NME. response is observed in the X'T' plane, and 
Consequently M (t) is the variable of interest. These equations 
are coupled^ and the nonlinearity is imposed through the time 
varying rf field envelope B(t), which appsa.rs as a coefficient 
in, the differential equations. Ihen B,j(t) =0, a situation 
which arises on the cessation of the excitation, the evolution 
of the magnetization M is governed hy linear differential 
equations* Therefore the BID signal may he thought of as the 
output of a linear system, whose initial conditions are 
related nonlinea,rly with B,j(t). This fact is brought out 
more clearly from the geoEietrical arguments to be presented ■ 

in section 3*^3* 

A solution of Bloch equat ions(3 *3) # in the presence of 
B,j(t) and field gradient niHfr''), is reviewed to sootioa 3«4» 
The problem is espacially relavsat in rf pulse tailoring for 
narrow band selective excitations. 

3*5 THE BASIC NMR IMAGIl^G- EQUATION (From geometrical arguments 
in the rotating frame of reference) ; 

A rd^oJTOUS solution of the NMR imaging problem lies 
with the Bloch equations, in spite of the model being based 
on macroscopic phenomena. However, it is possible to gain 
useful insight into the imaging equation, by considering the 
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geometrical argument of the process ing magnetization. The 
object of the present section is to show that the imaging 
equation may be arrived at from rotating frame arguments* 

The following assumptions are used in the derivation 
of. the equation ; 

t. The larmour equation w^ ^ valid during the 

pulse duration T » i.e*, S 3 ratem transient times are very small 

sr 

compare d to 

2* The pulse duration T^ << T.j 

The modulating envelope B|(t) is rectangular. 

ror the sake of simplicity, the case of a single direction 
gradient, G-^( z) in the z-direction is considered. The symbols 
to be used in the discussion are 

M^(z) ~ The equilibrium magnetization in the interval 
nz around i. 

?(x) ~ One dimensional spin density distribution 

~ J X (2:»7fZ) dy dz 

Y Z 

A -* : The cross sectional area, of the sample. 

Figure^ .2, illustrates the precession of the resultant magneti- 
zation. The precession cone is indicated by broken lines in the 
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f igure* Tlie cone degenerates to a circle at x = 0 whore 
The associated 8 (iUations are, 

L 

®res = ((^H )2 + ( 3 ^ 2)2 

tan & = /G-^x 

0 = T ^res ^ flip angle. 

At time T^, the magnetization has the orienta,tion 

determined hy 0 , ©■ and M^(x),the length of the magnetization. 

The X and Y components of the magnetization M and M respectively 

x y 

are proo®G'fcioiis. of M^(x). After the cessation of the rectangular 
pulse < the magnetization processes freely about the Z axis 
at a rate a w =^aH( x) . Ifence 

=:.M^(x) P^'(x) 

and = M^(x) " 

where P ^( x ) and P„f(x) are project ions of M^(x) at time T . 

X y o p 

writ ing Ef ( x) as 

M (x) = M^»(3:) - 3 My>(x), 

or p (x) = P^jKx) - 3 
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'the magnet izatio a M^(x) in turn may be expressed as 
M^{x) =Se(3t) 

■where K is a constant of proportionality relating spin density 
to magnetization# 

At the end of the pulse duration 

(3:) 

The subscripts indicate that the phasors M and P a,re in the 
x‘y* plane, and is implicitly assumed hereafter# 

The future evolution of the magnet ization is determined by the 
field gradient nH = G x. Thus, 

jyGxt 

JT ( t) = K P (x) ^ (x) e ^ 

The aggregate response from the whole sample is the magnetic 
mnment ^ 

, 'Jr: (t) = KA /f (x) fix) e ^ dx (3*5) 

X. 

The derivation neglects the effect of T 2 and T| on the assump" 
that 2 « T^ , T^ which is not valid for the longterm evolution 

PjL'x I d 

1 

of the magnetization* Its first order effect (when "(jJ- << ^ 2 ^ 
is, a multiplicative term ^ e f([cp. 
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Eberefore the response ^(i;) is given by 

w(t) = M ( /l’(z)f (x) e dx) e““^/^f(T^) (3*6) 

Ehis expression is cast in the form of the E^ourier transfer- 
mat ion, and is the basic UMR eOiuation. The term P( x) serves 
as an aperture function, -which modifies the effect of ?(x)* 

In the case of the rectangular pulse, 

^(x) = sin 6 cos Q (1-cos0) -j sin & sin 0. 

while the expression is fairly complex in the case of arbit- 
rary pulse excitations. A more detailed discussion is 
presented in appendix A. 

The i-m.aging equation suga-ests the following * 

1- The response m(t) is linearly related to<p(x), 
the system identification function, p 

2* The aperture function P(x). is nonlinearly related 
to the rf excitation ( t) • 

This is the nonlinearity referred to in the Bloch.' 
equations* As a result the solution of Bloch equations are 
necessitated to determine P(x) or the design of the pulse 
B^{t) to obtain a specific form for P(x)* This problem is 
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important in selective ezcitations, iwhere the resolution of 
the image is deterniined "by the aperture function. 

More generally, when the effect of the relaxation are 
spatially dependent, the ‘behaviour is no longer exponential, 
IThe resultant transverse magnetic moment m(t) is 

m(t) = EA J ^(x) f (x) L( x) e dx (3*7) 

^s 

where the function L(x), accounts for the effect of the 
relaxation phenomena. ' 

!Ehis equation may “be interpreted in terms of transfer 
functions, where (^(x) and l( x) define the system parameters 
to ,be identified, and m(t) may he treated as the impulse 
response of the transfer function E(x) ^ ( x) Ii( x) as shown 
in Figure 3.-3, 

The imaging equations may also he derived directly 
from the Bloch equations and is given hy Hinshaw [20] • 

The general. izat ion to three dimensions yields 

m(t) = E / 

O-s) 


f • 


1 - 
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3*4 A SOLUTION OP THB BLOCH EQUATIONS (perturbation teehnique 
due to Hoult) 

SectionS3«2 and 3.3 indicate tlie formulation of the 
Bloch equations in the rotating frame of reference, and the 
derivation of the basic imaging equation (3 ‘C) , identifies the 
aperture P(x) with the nonlinear component of the Bloch 
equations* 


Por the rectangular pulse^P(2;) is given as 

B 

, 71 / 2 ) 

ana S’U) = (, _ ,,03 C.9) 


for a 90° flip • The solution for arbitrary excitations may 
be obtained numerically using a piece— wi^ rectangular approxi- 
mation in the rotating frame (Gaprihan [t9l)« It is shown that 
gauss ian, sine and triangular modulations are good excitations, 
though not optimum* A more rigorous analjrtical solution is 
attempted by Silver et al.., [21]* However Hoult ’s perturbation 
solution [^7]» gi^ss an insight into the nature of the equation. 


The magnetization M(t) and M/t) a.re written as a series of 
perturbation terms of order yB-jCt), 


M(t) t= 

M»{t) = 

z z z 


(3*10) 
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bien as perturb at io a terms, and the Bloch 

j -,*£(?>)) = - aYB^(t)M/^~'‘^ 

(t)Im ('^*11) 

ith the Initial conditions =0 and = M^(r^), while 

he effect of is neglected* 

Phese equations are recursive, the higher order terms being 
determined by success ively lower order terms- The equations 
f#r the first few orders is given below, 

+ ih ^ 3YAH(r>)) =-*jYB.(t) M , 

dt ^2 

=> = 0, 

YB^(t) ImM^^^(t) 

"fc 

f YB^(t) ImM^^^(t) dt. 

-00 


and 


dt 


or 


M. 


.(2)i 


at ions yield 




dt 


(i 


dM^' 


,(n) 


dt 


!= YB. 


( 5 . 12 ) 



5 -^ 


It where /iiH(r^) = 0 , the magnot ization precesses 

about B|(t) or the X~axis. The flip angle for this case is 
given, by ^ As the effect of the perturbation term 

is in the product B 4 T (or in Bf for a given T ) the flip angle 

* * Jb^ 

0 may be treated as an indicator of the nonlinearity. 

The following observations may be made from the per- 
turbation equations : 

f 1 'l 

1. The equation governing IV ^ is linear, and is pre- 
dominant for smaller flip angles. 

. 2. At flip angles such as 90°, the transversal response 
is nonlinear, though the magnitude of the corresponding magne- 
tization is maximum. 

3* Operation in the linear region has the advantage 
of easy signal design as P( r) » where *3 denotes 

fourier transformation. 

4* the presence of a strong linear gradient, the 
magnetization at time T^ is negligible, and may be observed 

only with a spin echo pulse. This is due to quick dephasing 
by the gradients. 

5* Spin echo should be employed to observe the complete 
response in the selective excitation case, even if a. laoRge flip 
angle is employed, as the primary response is linear. 
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6 * !I3ie excitation should not possess sha.rp edges- 
Hfhis increases side lohe power la the primary response^ 
which reduces selectivity* 

7 * Higher order teirnis of F( x) ( /v M( x*3!p) ) hear no 
resemhlance with the spectrum of the pulse* 

The latter ohservatioog are based on the assumption that th© 
linear response is significant. Hoult [ 17 I indicates a 

ratio for a rectangular pulse 

2 

^ 

“8 

which for 0 = m /2 is 

A comparision of the Bloch equations with the geome- 
trical interpretation is given in appendix A* 

Practical considerations for signal design are given 
by Carpihan [tg]* Selection of an excitation is based on .the 
behaviour of PCz)- The various parameters describing the 
figure of merit on l‘(x), for selective excitations is shown 
in Figure "(2*62} • 
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3 .5 SYSTEM EQUATIONS FOR DIFFERENT ELSMSNTMY EXCITATIONS ; 


The generalized imaging equation (3*8 ), describes 
the response of the MR system to a single rf pulse* The 
function F(r^) may be obtained as a solution to the Bloch 
equations. The present sections formulates the equation for 
different basic excitations. They are, 

1 . The rf Pulses ; 

a. 90° and 180° broadband 

b. 90° and t80° narrowband 

2. Gradients 

a. Readout gradient 

b. Gradient pulse 

C. Gradient reversal. 


all of which are treated in sections 2.6 and 2.7. Whereas the 
rf pulses, affect the response M(t) through the aperture 


function F(r^) , the gradients modify the response through the 
phase term e ^ ^ * . It may be observed that the image 

resolution is affected by both functions, though primarily 
by the gradient* The equations of the excitations follow: 

1 »a 90° and 180° Broadband Pulses ; 


This corresponds to F(r^) being a constant over the 
sample space* Geometrically, this implies,B^ >>G^, T^«T 2 a 
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!E| and 1/ 5 ‘(r^) = -3, resulting in 

s(t) =1: ^ ^ (r^)e"'^^'^^ Q3Y*3-(r ) .r ( 3 - 15 ) 

7 

1 .b 90® and 180® narrow band pulses ; 

The pulses are primarily designed for slice selection. 

ITbe general condition on P(r^) is shown in Figure 3 * 4 , This 
corresponds to 

B^ C5i.G2.2: for z < I, T^ « T^ ,T2 and 

ISC 

and a proper choice of B^(t). For rectangular pulses, F(z) 

(for a single dlroetioi^^ shown in section 3*2, and the requisite 
pulse duration is given by the equation, 

1 

(y(B^^+(G^z)^)^ Tp = ^or a n pulse (^-14) 

2 *a Readout Gradient ; 

This is applied foll.owing a non -selective x£ pulse, to 
effect the readout of the FID* The imaging equation is 

m{is) K f f (r^)F(r^)l(?^) 


dV (3.15) 
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similar to equation (3»8)- It ig essentially a free eirolution of 
th.e magnetic moment, in the absence of an excitation. 

2.b Gradient Pulse : 

IThe duration of the pulse is shorter tha.n t/yG^^ so 
that the magnetic moment is constant o-ver the interval t and 

m(t) s= K f ^ (t^) ^ dT. (3*16) 

^3 

fhis is the Issy equation for Courier Zeugnsitography* 


2 *6 Gradient reversal » 


Broadband spin echo may be obtained using these pulses. 
Ihe gradient G(x) ( in the x direction), is reversed at tline 
after the application of the 90*^ pulse, fha evolution of 
the pulse is given as 


m(Ijj) = M / f (x) 


e ^ ® dx 


aad m(t) > KA/f(^)o 


dx 


/ ;JYG_x(2!E-,-t) 

= My ^ (x ) e ^ ® d 

for t > 

at t »(2fjj) = KkJ" f (x)dx., 

corresponding to a rephasing. 


(’^. 17 ) 



IDhe equations for different modalities are obtained 
as combination of these excitations* Appendix B lists the 
equations for some salient modalities and for the recovery 
techniques namely the SR and IR. 

FIB iin) RWR ESCEI7RR : 

The precessing magnetic moment m(t), described hy the 
equations (3,7) and (3»8) , 5s picked up hy a coll- The 

induced emf is given hy Faraday’s law, expressed in the form, 

= ( 5 . 18 ). 

using the principle of reciprocity, (Hoult [22])* 

§1 denotes the magnetic field produced hy a unit . 
current in the coll, and is assumed to he homogeneous, over 
the sample volume, and ffi^(t) is the magnetic moment in the 
iahoratohy frame of reference* It is related to m( t) hy 

jw.t 

a(t) ^ m(t) e ® . (3rl9) 

T^ing the equations (3*i8) and (3.19)# l^he emf is written 

as , , , ' - 

;jw_t 

= - j B^w^m(t) e 

for, the Bandwidth of m(t) is small compared to (the narrow 
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band assumption). The coil is linearly polarized and the 
FID le given by 

t 

s(t) = Re (-3B|W^m(t) e , ). 

"b 

= B^w^Re(m(t) e ^ )* C^*2o) 

Figure 7*6 * indicates a typical 1-D image distribution. 

OJhe image is assymmetric about z=o* as^ Cr) will in general 

be arbitrary, while F( x) is complex. Consequently the 

demodulation the FID is done at a frequency reference w 

different from w„, 

o' 

"w 1 w - Y ( ^J.r^ ) , 

max 

i..e. the upper side band. 

In actual ME. systems, the irradiating frequency is kept 'the 
same as the detector reference frequency w» The rotating ■ 
frame analysis is then done at an angular frequency w, and the 
residual field along the Z-^xis is given by 

H(x) = ^2-2: > (Wjj - -w). 

5(w - w)t 

The equations(3 ,6-5.8) are modified by e . The 

magnetic moment in the laboratory frame of reference however. 
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Is the same as in the previous analysis. This is because 

a change to the stationary frame is effected by a factor 
i-wt 

, instead of e ® . 

The equation ('5 .20) , is rewritten as 

/ 3 -.-i-Tr/o 

s(t) * A B|W^Ee [{ y M(x)e ^ dx) e e ] 

"" ( 3 . 21 ) 

where M(x) = P^y.(x) M^( z)l(x,t) is in general complex 

M(x)- i Mg_(x) - 3 Mj(z) . 

Evaluation of the real part yields 


s(t) * A w^ ij Mg^(x) cos ((w^+Y “ V2) dx 

Mj{x) sin ( ( w^+yG^x) t - Tt/2) dx] 

Quadrature detection with reference signals cos wt and sin wtj 
( Including low pass filtering) gives, 

B^{t) = — [ ^Mj^(z) cos|^ (w-(w^+Y&2.^))'t+V2j- dx 


x) 


sin I (w 


« (w^+Y0^x))t 
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AB* w 

and Sg(t) = 



z) (■w-ftf^+yG^z) ) t-!-ix/2 dz 

+ Mj(z) cos ^ (w~( w^ -!- YG^z))fc+ J dz ] • 


Alternately , 


s^(t) = 


A B^w_ ^n-rr/O 0 (^^«+YG'^X-'w) t 

--^Re r/M(z) e 3^2 e*'' ^ dz] 


ij M( z) 


and Sg(t) = la M(z) 


3 (w +yCt z-w)t 
e ° ^ dz] 


(^. 22 ) 

Considering that in an ezperiment the detector reference w 
3s fized,(w^ -w) ^ 2 ' seen that s(t) 

jsj't) may he Sourier transformed to obtain M(z). The phase 
introduced is in genersl different from V 2 , due to the 
contribution from the receiver and has to be corrected. The 
conceptual schematic of the receiver is shown in figure 3 -5 • 


A special cs.se of the equation (3 .22) , arises when 

F(z) (same as F ( z) )= — j, iCx) constant over z and w = w . 

Xj O 

Here s^(t) = - Re (M^(z))] 

and Sg(t) = - \ (Mq(x))]. 

This is illustrated in figure 3 * 6 , and it corresponds to a 1-1) 
projection of the spin image f (z,y,z). 
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The other situation encountered is PCs:) = F^( 3c) - y) 
and X(x) constant over x. This represents a selective pulse 
excitation, where ~( x+x^ )) = - conjugate 

skew symmetry), while ^^( 2 :) is constant over P( x+x^ ) . x^ is 
related to w (the operating frequency of the transceiver) by 
■w = • 2!he quadrature component: it' C^'n 

Sg(t) = Re [IJ (PCX))] 

In most modalities the signal is observed in the presence of 
a rea"'.out gradient, describing a combination of the two cases 
considered above# 

In addition to the PID, noise gets coupled to the 
system at various stages of the receiver. Its effect degrades 
the quality of the spin image obtained, and has to Hiinimized. 

In most HMR experiments signal to noise ratio is sufficient 
to deteraiine the sensitivity, while in some imaging applications, 
noise statistics is warranted for developing optimum receivers, 
and for image enhancement. 

Figure 3.6# shows the noise power spectruni along 
with typical a signal spectrum. The source of noise is pre- 
dom.inantly thermal, and is assumed to be white Gaussian, 
around the frequency of ope rat ion. Considering the receiver 
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l^andwidth. to be larger than the signal bandwidth, the noise 
is assumed to be a bandpass, white G-aussian process* Denoting 
and Jig(t) as the quadrature noise components, 

(the cross correlation) is 0, as the spectrum of n(t) is symmetric 
about the carrier. The two noise processes are thus, uncorre- 
lated, and can be treated as statistically independent white 
Gaussian, processes with variance 

The fourier transformations 
00 

= J n^(t) e”^"^ dt 
—00 

00 

ir^(w) = ^ ng(t) dt 

—03 

preserves the correlation and the Gauss itiQi property. 

¥ 

For a real function M (w(x)), s„(t)) is real and 

o ^ o s 

the reconstructed image is 

\ (w) = M^(w) N(w) 

where K(w) « Re (H-(w)) -i- 110(11 (w)) 

0 s 
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In general, where M^{w) is complex, the Image may he 
defined as M^(w) and the noise will be a Rayleigh distributed 
random variable for any w. 

The SFR for the experiment is evaluated, for the 
corresponding bandpass signal f(t) * 


The signal is assumed to be time limited (by observation) 


Por M( x) = ^ constant , 


I ./2 


f(t) = A Re [ j (e ' ^ dx) e ° ] 




-L/2 


G L 

sin Y 


3— t 


= w B^ A L 
o 1 


G L 


cos w^t (t. 23 ) 


Y 


The signal is gated 'on’ from t = 0 (actually T^, which is 
small compared to 1 /Y^^p ^2^ * 

The energy B in the envelope is, 
s 



asp o 

/ V at 

0 (Y -fc )2 
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2 2 2 

= Si V ^ ^2 


/m * ^ for strong 

^ ^2 • yOrJu 


gradients) . 


(!!. 24 ) 


The maximum SRE. obtainable then, is =v 

where is the white noise power spectral densitjr^if the 

signal is known apriori^ The ratio can then be used as a 

representative figure of merit, as it is based on the 

assumption M(x) = general the spin density variations, 

and the T^ , T 2 variations in the sample are small. lEhough 

the variations ( 9 and T 2 (?^)) carry the 

imaging information (as contrast), the bulk: sensitivity of 

the system is indicated by E_/E . Additionally, this compares 

s o 

well with the derived, for a pulsed HMR experiment, by 

Barast and Anderson [ 1 1] • 


To sum up the system identification problem is, ela- 


borated in figure 3.6. 


CElLSTm 4 


SISTIM DSSICtIT cons idbrat ions 

4.1 SYSTEM SCHEMATIC 

A system identification problem was defined for HMR 
Imaging in the previous chapter. Figure 4,1, Indicates the 
conceptual schematic of the system. The present chapter 
considers various aspects of the system, with the view of 
developing a design strategy. The block schematic of the 
system, is broadly class if led into five categories, in keeping 
with the conceptual schema. They are 

1 . Sample coupling 

2 . rf transceiver 

3. Signal aquisition and averaging 
4* Magnetic field dontrol 

5 . Sequence control 

Saiiple coupling, is the centre of the system, acting 
as an inteTface between the sample under study and the 
outside. This consists of 

1. Main field coil system 

2 . Gradient coil system 

3. The rf coil. 

There are briefly described below* 

1 .■ MAIN FIEID con SYSTEM ; This is an important part of 
the system. In as much as, the magnitude of the field H^, 
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(161:63:111131.63 ’ulie operating f x'eq.u8ncy The spatial and 

temporal fluctuations over the sample, tantamounts to 
degradation of system perfo3m3ance. In high resolution KMR 
spect 3 ?oscopy^ shoi*t tern tempora,! fluctuations, and spatial 
fluctuations are maintained within 1 ppm of the main field 
(which is typically 15”25 EG for protons and higher for 
other nucleid'S-s) ♦ This constraint on imaging Systems^ is 
less Stringent., due to a deliberate line hipadening introduced 
hy the field gradient However the stability requirements 
on the power supply remain owing to ths large sample volumes 
involved* "Where imaging times are large, long term field 
drifts are traced using a known sample and frequency locking* 
Hoult indicates a design of an electromagnet for medium 
fields (^1*4 K:g), using Helmholtz pairs on a sphere [23]* ih© 
homogeneity is reported to be 3 ppn over a sample volume 
of radius tO cai* 

2* GHAHIBHT COIL STSTSM 2 Three pairs of gradient coils 
.are used for and Z directions respectively* The first two 
are generated using a saddle shaped coil geometry shown in 
Figure 4..2. This is known to give acceptable higher order 
spatial field nonlinearity* The Z direction employs a Helmholtz 
pair.., co-axial with the main field. Typical gradients 
acheivable in large magnets are,-^- 0*1 “0*5 G/cm* For saddle 
shaped coils higher magnitudes are accompanied by greater 




FIGURE ^ i BLOCK schematic OP NMR IMAGIWG SYSTEM 
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nonlinearity or for a higher coil separation, higher currents. 

As the gradients are generally switched, the associated transients 
are large for high currents. Thus a compromise is sought "between 
current, gradient magnitude, switching speed and nonlinearity. 
Switching times for a gradient of 0.2 G-/cm (at 30A current in 
whole "body imaging ) is 1 m sec. 

3 * yf. 00 XL This couples the sample with the rf excitation. 

A single coil geometiy ' is employed, usually of a solenoidal or 
a saddle shape , though single turn loop pairs are known to 
be used, in crossed coil arrangements. As these coils produce 
linearly polarized fields, only half the current is effective 
fCr sample excitations. Tie performance of the coil is measured 
in terms of rf field uniformity, and couplins efficiency 
(related to SIR). As these requirements conflict, a tradejoff 
is employed .Tor Cx£m\plej,a.:roleno id has a lOdB SNR advantage 
over a Saddle Shape of same dimonsions, though the field uni- 
formity is poorer. J‘or..a radius ratio (sol. to saddle shape) ^2, 
the SHE advantage is still ^ ?dB, for the solenoid. 

The following section presents the design considerations 
for the sample coupling, rf transceiver, signal aquisition 
and averaging subsystems, while magnetic field control and 
sequence control are not within the scope of this thes is . 
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4*2 GElJERiO:, DESIGH COiTSIDSRA'TIOL'S ; 

The applications of MR to clinical (e.g. proton imaging 
of organs) and chemical ( r n-^'lvo metabo3.ic dtudies using 
Contexts are many. These areas call for different design 
considerations, owing to the complex Interrelationship "between 
the observables (P(r^, T^(r^), ”^^3 imaging 

environment. For exam pie, detectability of lesions is more 
dependent on T^ [ 6 ], while imaging of various organs are 
achieved b3r mapping^ "(r^) only. The modalities in these two 
situations are different, the former req.uiring methods such 
as inversion recovery* An overall design perspective is 
called for to sort out - the specificities and common points, 
amongst different contexts. 

Kaufman and Crooks [7 ] indicate general corisiderat jons 
in imaging, from the point of view of image quality, SHR and 
imaging time. The choice of field strength H^^Cor the operating 
frequency w^^) 8.nd modalities f or at^lering different goals are 
discussed. Wagner and Brown [24] propose the problem as an 
unified approach for all medical imaging methods* The 
"detection (i*e- image aquisition from the transducer to 
digitization), and the display stages are distinguished. 

\/hile the former determines the best possible image quality, 
the latter interfaces the data with a human observer. The 
analysis of the detection stage involves the concept of an 
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ideal observer, defined in the statistical senoe. Image 
quality is related to SUR in terns of a contrast~detail 
curve of the form C^cl?^Q=Constant , for the X-ray CT* C, d and 
Q denote threshold contrast, lesion dostmeter and SRR (or 
Noise Biuivalent Quanta, in X-ray -context) respectively* 

Ihis approach sets bounds on the detection and display stage 
performances* 

Ihe present discussion assumes, Image quality and 
^jiiaging tome as ’object ' parameters, for system design, along 
t^he lines of Wagner and Brown [24] • These are translated to 
system specifications (for the detection stage - as the 
discussion is on imaging hardwai-e) , namely SNR, magnetic' 
f ield Field Gradient Transmitter pulsed power W^ 

and the modality* Two specif ic at ions are also determined by 
second order effects such as^ 

1 • Instability and Inhomogeniety of 
2* Nonlinearity of gra,dients 
■ 5 * Non-uniformity of rf field* 
over the sample. These are implementation specific and are 
not discussed here* 

Image quality, is defined in terms of, 

.1* Resolution s6^ (The smallest distinguishable section 

in an image) . 
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2* Contrast : (noimalized relative int 9 i]eity of 

o 

tlie object I'rith respect to back- 
ground) 

Texbure is also a parameter, and is seen to depend on the 
jm'aglng hardware nonunif oimities [7 ]• However this is 

not discussed here* 

Resolution itself is defined in terns of, 

1. A point pair 

2* Contrast ( 4n noisy images) 

5 » temporal, conterfes* 

Ifhis parameter sets a lower bound on SHE, and is important,- 
when signals are very weak and temporally fluctuating (as in 
whole body scanners and in WLR microscopy [ 25 ]) * The depen- 
dence of the object criteria on the detection stays is quanti- 
fied as shown in Table 4*1 , while the general design conside- 
rations are summarized in Figure 4*3* This shovs ^ that the 
specifications of the sample coupling, transceiver and signal 
processor subsystems are affected by the dimension of the 
sample, image quality (resolution and contrast) and imaging 
time* These specifications are discussed in the nezt section. 

4»3 SYSTEM SPEC IF ICATIOHS ; 

The first order design equations and the assumptions 
involved therein are given in Table 4*2 . A comparision between the 
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two systems (a multi-phase bio-system jmager and a clinical 
imager) is presented in Table 4*^, in the form of a case study. 
The response of a single 90° pulse is considered. The sensi- 
tivity of the imaging experiment is determined by the signal 
to noise ratio. This specification is discussed to start with. 


1 4 SICTAl TO NOISB HATIO ; 


The equation( 4*^ )(Table4.2) shows the dependence of SNR 
on the various system variables. This is defined as a ratio 
of. signal energy to noise spectral density, as shown in section 3*6 • 
Tlie SNR for an equivalent pulse resonance experiment is given 
by equation.^ 


SNR .== K. 




7/2 


A 

2 




G-, 


SA 


These equations are valid at frequencies^ where inductive 
losses are not predominant. An estimate of the relative loss 
is. quantif ied by the figure where R^ is the resistance 

of .the RP coil, while R^^^ is the equivalent series resistance 
accounting for the dissipation in the lossy sample. An expre- 
ssion for the induct-iv© loss is gr;en in equations (4*4 and 4 . 5 ) 
for solenoid and saddle shaped geometries respectively, Hoult 
[ 26 ] ♦ Consequently the signal to noise ratio increases with 
a power of 7/2 for w^ « w = R^) , (the frequency at which 

* all equations are listed in Table 4.2 
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= Eq) t and a power of 2 ', for w^ » w = S^) . Por 
example, the case study (Table 4 . 3 ) indicates 


f(®m = ^c^ = ^ =3/2at f^ = 15 MHz 

K 

f{\ = Rp) ^ 0.89 MHz and ~ 17 at f^ = 6 MHz 

respectively. These values are computed using a lossy sample, 
of loo mM KaCl with a resistivity of 0.85 ohm-m. The other 
factors increasing SHE are the sample volume, V signal 
averaging gain G-g^» sample characteristics Kj and coil factor 
Kcoii* Solid or solid-like samples with lower relaxation time 
T 2 » have poorer SNR compared to liquid samples. A detailed 
discussion on Gg^ from the point of view of modalities, is given 
in in Branner and Ernst [54]» 


2. MAGBNTIC PIBlD H^ : The choice of field magnitude depends on 
a* Availability of magnets for the task, 
bi Power requirements on the magnet power supply, 

c. Complexity of electronics at the operating 
frequency 

d. The cost of performing shimming and stabilization. 

These govern the economics, and tethnology of choice, while 

the advantage of increasing the field is based on the relation 
7/2 

SNR a (Hq) . This is not valid for high field magntidues and 
the detering factors aire 
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a* Self resoaaace effect of the coil at high freiaenciss . 
Ihis frequency decreases with increasiag coil dimensions, owing 
to distributed capacitances Cd* There are also associated 
losses in the sample dielectric proportional to Cd. 

Por example, in head imaging Cd = 20 pf, (Cd 2pf/cm 

diameter in sol&ndids), which for an inductance of l=1viH 
leads - 25 MHz or = 6 KG. The dielectric losses 

may be represented as an effective resistance R^, where R^/R = 

T Q , and T, Q, C^^ and C are loss factor, coil Q, distributed 
and tuning cajHiGitances respectively (Hoult [26]). These may 
be minimized by reducing Q (or w^) or by increasing the self 
resonance frequency by choosing coil of' smallest inductance ■■ 
for given dimensions with acceptable homogeneity* 

b. Inductive losses in the sample are given by 

'2 2 , 
m«- ^ b'^ 

R^ Kj j 3' unifom spherical sample of 

s 

radius b and resistivity 9a* expression is derived assuming 

negligible skin effect in the sample. Por a lOCa^, Nad ■ 
solution, with b = 15 cm, = 0.03 ohms at 4 MHz, 0.5 ohms at 

15 MEz.and 3 ohms at 40 MHz. At higher frequencies, the 

' 3/2 

inductive losses are predominant, and SNR or ’(w^) . 

c. Radiative losses in the coil may become significant 
for wavelength >'^1 and has to be reckoned with in whole body 
imaging, where Ii'^6m^ wavelength ^ 2Caa. 
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d* T:j is kao'wn to Increase with implying a larger 
imaging time using SR and IR techniques, discussed in chapter 2» 

Consequently the field cannot he increased indefinitely* 
iPypically the operating frequency is deteimined hy the sample 
losses, and hardware costs. At frequencies w^ ^ w (\=\)» 
the SRR improvement is marginal while,- the technological 
factors may hecome overwhelming. As an example, superconducting 
magnets for medium fiel^ 6EG) do not give significant per- 
form^ce improvement* locher [15], indicates the use of water 
cooled magnets at H^=!l.4KG, for conmiercial applications. In 
Tahie 4«5| the fielcisare chosen at* 5*5 EG and 1*4 EG for the 
two Systems respectively. A value of .4EG is conservative, 
due to the small f g used. Realistic operating f ield range of 
5“4iS- is suggested for medical imaging hy Eaufman and CroolJS [ 9 ], 
by assessing the relative SER contrihut iocs of various organs 
of the hody. 

^5..■ Pmi) GRADIENT AED RBSOIUTIOE ; 

, The resolution of an image increases with the field 
gradient as shown hy the equation 

w = Y (G{r^) • Ar^) 

However Increasing the gradient heyond a limit has the dis- 
advantages that follow, 
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a«, SUE per unit freq.uenc 7 decreases due to signal 
energy spread over a largerbandwidth, for given noise spectral 
density. 

i 

"b. 12 '^ ~Y &1 '^©Greases, rendering signal aq.uisition . ^ 
difficult, as attenuation "bet-ween adjacent samples is e 
This requires signal averaging, thus increasing imaging time. 

This is disadvantageous, where sample movement is unavoidable. 

Table 4-45, indicates a calculation for resolution. The 
intrinsic frequency resolution limit considered is i^-V^I/T , 

■; W 

where T^ is the observation time window for the PID. In imaging 
experiments T^ < T 2 . An. optimistic estimate of the smallest 
volume observable (for a 0 dB) , gives 

V ^ 3x10"''^° m^. 
and V '^ 7 . 8 x 1 o“^m^y 

for the two systems respectively, with T^ = 500p See or = 2 KH 2 « 
It may be noted that E /M is the same as the ratio of signal . 
power in the volume to noise power in bandwidth av • This 
compares well with = 7 > 43 xl 0 ‘'^° m^ by locher [I 3 ], obtained 
by extrapolating SHE from the tables therein* Assuming a 
cube of volume ( ^Y) = (ax)^, the linear resolution for the 
two systems are 

n X r-o 0*7 mm at SM= 0 dB 

and X ^ 4.3 mm at SKIfe 0 dB respectively, while in 

generai ax oc SHE* That is at SEE = 10* Ax = 7 no®. 



W* House [ 25 ] prosents similar calculations to show 
resolution in microscopy context. At 600MHz, a cylindrical 
pill "box of length 20va3t and radius 2 ijm , yields an SHE* of 10 . 

A frequency resolution of iHz is assum.ed for calculation. 

A more realistic estimate of the resolution, has to he 
defined in clinical imaging context, where contrast and 
resolution are related, as in detection of lesions in an image 
■field. IPhe perfoimance of the detection syst^mi has to he 
evaluated using a conlniast -detail curve of the hind derived hy 
Wagner and Brown [24]- This is also, appropriate in imaging 
other nude ide distributions {e»g., where the inherent 

sensitivity in terns of is poorer. 

4,. BIGHAl lETBl AND TRANSMITTED T0?JSR ; 

These factors are based on the dimension of the sample 
in as much as the coil dimensions depend on it* As ouch, they 
are dependent variables^ specifying the receiver gain and 
power amplifier requirements. The dependence of hardware 
specifications on this is shown in Table 4.4* The values given 

in Table 4 •'5, are representative* They are, 

7p7/turn and 0.84mV/turn 

250 - 1 OOW and 200W-1 .45hW 

pulsed ms pulsed rms 

for T 5-10 for T 500-100 

pSec.P pSec.^ 


Signal level ; 
at t=0' ■ 


¥. 


P 


* POWER SNR C^ 'oAb) 
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for the two S 3 rstems respectively* As the numher of turns in 
the coil (n) is specified hy the coupling circuit configuration, 
the emf induced is expressed in volts per turn, fhe transmitted 
power (pulsed rms value) is derived by assuming that the power 
coupled to the rf circuit is dissipated wholly in the sample 
coil and the sample itself. Squat ion( 4 ,^ indicates a 1 /T ^ 
dependence between pulse duration and and a 0 ^ dependence 
between the flip angle and ¥ » Consequently, it is unrealistic 
to produce broadband 180 ^ pulses in high gradient experiments. 
Spin echo, is usually effected by the' use of gradient reversal 
(which again is speed limited due to switching rise time) * For 
example, ¥p for a 180°, 50jiSec rectangular pulse is 6 k¥, while 
for a 25 pSec pulse would be 24h¥. However^ as 180° pulses are 
used, to produce spin echo, its duration may be greater than 
50hSec* fhe power requirement ooatin,n.es to be large , consi- 
dering dissipation in other parts of the transmitter circuit. 

5 .- SIONAI AOQtriSITIOH AHD SAmUJg V 

The FID after demodulation is recorded digitally. The 

minimum Sampling rate F„ required is deteim.ined by the Uyquist 

s 

criterion, F >2 (B.¥. of the signal), while the number of 

* s 

samples per unit bandwidth is determined by resolution consi- 
derations. The intrinsic resolution is assumed to be 1 /!E,» „ 

-tf t/ 

where T^ is the aquisition window. The FID is, ® 
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(actually it is more realistic to approximato ae 

"^g = ^inc ( ^ ) and = y"^L^ * ®ils corresponds 

' X 

to the signal for a 1-D projection for a gradient in one 
direct io.n« 

2 ^ 

Calculation of the time 1 at Tfhich ktl > (the 

xn® noise voltage), yields a lower limit on the observation 
time^ Table 4*3 » indicates a sampling rate, sampling interval 
and the aquisition time T(Vg = <n^> ^ of 70 IHz, 15pSec and 
600pSec for a 50 pSec decay signal and 400EHz, 2.5pSec and 
tOOpSec for the two systems respeetively* Tho contribution 
of the receiver noise is not accounted for and hence this 
gives an optimistic estimate. 

This indicates a total of 40 samples in the observation 

duration for both cases. In reality the number of samples 

efc#?<lutirod in a single experiment is more, i.e., the observation 

time may be more than 2, defined above* This is because, 

2 - 1/2 

though the signal level goes helow the noise level <n > ■' , 
tW total SIR of the experiment (or SIR over the frequency 
spectrum) is within acceptable limits. The evaluation is 
shown in the preceding section* 

Consequently the aquisition time T^ti» limited by if 
pulse repetition rate (or modality e.g*, IR and SR sequences). 
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a? 7 Picall 7 for = 50mSec, may be as large as 25mSec, 
as shown, by Locher ['^3]*, The minimum number of samples required 
is determined by the ratio, K = yG-^/2nA'i> or E = ^ /nsz (where 
^ resolution limit in frequency and I, the length 
of the sample). For the two systems under study K (maximum 
required for a resolution at SHR = OdS) is 30 and 60 respectively. 
More number of samples or pixels in an image for a given T^, 
does not improve this resolution, though may be convenient for 
viewing and interpretation. 

4»t4 , AN SYSM FOR MOITIIS^ 

The first order calculation, based on the expressions 
given in Table 2, are perf oimed for an imaging system used in 
multiphase bio-sample study. The results of the calculations 
are summarized in the first column of Table4*3. The assumptions 
involved are listed below • 

1* The imaging requirement given in the beginning of 
Table 4*3, cover major bio-samples of interest. 

2. A spherical sam.ple of resistivity 9 ^ = 0.15 ohm - m 
(lOOmM NaOl solution (Hoult [26])) is used for evaluating 
inductive 3.osses. This yields a value of on the conservative 
side • 

-,5 3 

3 . The effective volume of the coil (<-^5x10 m"^) is ten 

times the total sample volume, ensuring good homogeneity in 
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the saddle shaped (BerJaardo [ 29 ]) and solenoid coils . 


4*- 2!he iron magnet available, generates a 15KG field, 
at a current rsJ 3 O- 4 OA, while the saturation current is less 
than 50 A * . 

5* hnaging time is not a premium, as the sample movement 
is restricted* 


The design values calculated namely, 


H. 


3.5 m 


Vg(t== 0 *^) = 7 pT/tum 

¥ (180*^) == 300 ¥ {dissipation in the coil and sample only) 

ir 


W. 


avg 


1 . 0 ¥ 


and Tg 15iiSec. 


are ImpleEientable with existing pulsed MR facilities 


The gradient required for 1mm resolution (corresponding to 
''j 2EHz) is *'>' 5G-/cm while with a 200 Hz resolution is 
0*5 G-/c5m* The hand spreading due to a gradient of 0.5G-/cm:' 
over 5 cm of sample is lOEHz or the decay time constant is 
20 piSec (for an equivalent exponential)* The limiting factor 
on resolution is then the signal to noise ratio* The previous 
section Jndicatss an SHR of 0<© at 1 mm resolution* The 
performance may be riiproved by signal averaging at the cost of 
imaging time. Repeating the experiment ten tim® gives an 
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SUE. of 1 OdB af 1 mm resolutioi’i* Th.e ampere turns required for 
the gradient (0*5&/cm) is 1215 for a Helmholtz pair of R/L=1 

■ at 

and t9A5 for H/L=l/2* The ratio 2R/L=1 is optimum in the sens® 
that, the second and third derivatives at the centre of the 
co-ilacpo equal to zero. 

The signal to noise ratio, = TOdB for Tg* = 2.273:10'”'^ 

Sec and 55® for Tg-^OliSec, for the vrhole sample excitatioi^ * 

This corresponds to the PID from a single image project ion - 
A set of.sfe pro jection in a plane are tractable for recons* 
traction of the image, considering the soft-ware comple 3 ;ity. 

. The resulting is 77® (taking into account signal averaging) 

SO 

The corresponding (0® SWR) resolution is AY=5,.32xl0 m-"^ 

( or n^. 3 >=lmm) • In practice SRR < iioise and signal 

, handwidths are not matched. This may correspond to an SHR • 
degradation of upto lO® (after accoimtii^ for receiver noise). 
(Cohsidering the complexity of selective excitation and 
gradient s-witching, projection reconstruction is a plausible 
mo dal ity f or the appl Icat ion) • 

^ plSURE 4.1 ,p-7d. 
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4.1 

imm QUAIiITI MH) SYSTEM PARAlfflTlRS 


^ nnwfPR a qfTi 

SPATIAl RESOLUTION TEMPORAL RESOLUTION UOIliKAbi 


G-radient 


T 


2 


SWHM OP 
P 


SNR in 
Noise 
Limited 
case 


1 .Sample Motion e.g. P f-r'^ 
jfandom t'luctuations 

in sample T^ 

T‘ 

2 . Drift in v 

detection 

System 

Contrast 
Indue ing 
agents 


IMAGING TIMS 


Modality 

Signal Avg. 
reciuired(-f ) 

Computational 
power ( -) 

Saaaple 

movement 

(-) 



89 


TABLE 4.2 


FIRST ORDER DESIGN EQUATIONS 


EQUATION 


REP. 


SNR : 


7 /2 

'“o = ^*2 


L^.tL^a 
^0 


sA ^ ^ 


‘king equation, 

T. 2 2/A.TT^2^ * T3 2 

■H _ :2jk.'^Lj?2_ !l. 

o “ T^ * R 


Hoult 

[26], 


Imat 

[111 


( 4.2 ) 


ASSUirPTION 

a. is constant over tLe 
sample volume 

b. Self resonance, radiative 
losses in tlie coil neg- 
lected. 

c. Effect of finite Tj^ (the 
repetition interval is 
neglected) 


re 


= sample losses neglected) 

= R^+R« ( inductive losses 
° “ present ) 


general 

2 p 2 

^ (a Wo"^^}- P 


2 

0 


« a 

o 


1/2 


=> negligible 
losses 





>> 


a 

o 


1/2 


losses 


predominant 


SAMPLE LOSSES 

°d 

Lectric Re = tQ R 

EtW 2u 2i35 
J • Tl i 0 "^O 

active Em = —st— ^ 

^s 

a Effect 

o 

(4.5) 


2 P 
a u 'W' 

o 


Hoult 
[26] . 


a. A13- dielectric loss 

effects are luaiped into v 

h. Electricl field due to 
only 

c . Sample modelled as a sp- 
here of radius CT:^nd uni- 
form resistivityCp^. 

d. skin effect neglected 

e. Losses are lumped to an 
equivalent resistance Re 
or Em. 



equation 


ESP. 


ASSUMION 


. 00 IL PIQUES OP IffiEIT 

B, 2/R^ 

. COIL GB0M5TRIBS 


a. Solenoid 


. n 

|®l|o = - 




B _ ZSfiCS. 

V 3 I 3 


^UU ^ 

0 


^ = nW:i 


1+0.9(Ey 


R and L^are in cm. 

2 1^0^^ 

/^c = ia^TR^ 


n>>t 


Hoult 

fee ] 

Henney 

fe? 1 


1 . n the numher of turns 
o 

is large for a solenoid. 

2. O' is a proximity factor, 
rela.tlng to nsighhouring 
field effects 

3 . distributed effects ( c^) 
are not considered 

4- }iy = 1 

5 . The width of the saddle 
section, (n-‘1)3r < L/IO 


JUI 

for and coil separation 

3r. 


R«, = 


v^2 2 2, 5 

0 ^o 


30 Vg(R +(1^/2)'") 


(4.4) 


• Saddle Shei-oed 


1 3 


Bi I = (s^+s^)sin^ Ginzherg 

^ ° 4R [28] 

l-f(5f)2 


■ crn?»( 40R+4L ) 60n ^ 


I- = ( 2 B^) pi! 0 

,2 2 c 


rl 1 

, 2.o"2x2 


I^c ■ 20^ W2 " (V3 

L 

"or ,(n-1 )3r< |g- 


+ 1 ) 


md 0 = 120 ' 


( 4 . 5 ) 


cont inuted 
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EQUATION 


EBF 


ASSOMIOITS 


. TRANSMITTER POVSR RBQUIREESNTS 




B 


Hoult 

[ 22 ] 


1 


for a 180° pulse 
nd ¥ = ¥ 

'"avg ""p Tj^ 


(4»6 ) 


. PULSE DBRING-INO TIMS 

important in single coil rf proles) 


o 


. PULSE LURATIOl 


(4.7 ) 


roadland 


T = 
P 


T(}^ 


( 4.8) 


arrowband rectangular pulse ; 
inc(Y(Bf-i-(U^^x)^)^Tp=^|-^ 4.9) 


. INLUCBD eraf 


=s(t) 


t=o 


=w £^Y 

o I ■ 0 


(4.10) 


i. All power is dissipated 
into the coil resistance 
2- Only balf the current in 
the coil goes towards 
gene renting the rotating 
f ield. 


a. /^2 


. ORALIBNT 
i>f=4257.6 Hz /cm 
if= Af 6x0 Hz 


0 0 0 0 

*7 

-or resolution : ( 6z) ^=( 4. V) 
u I(27(R) 

=2 71; L=R Saddle shapped coil 

■^lO L=2R Helmoholtz pair 

(4.11) 


a. An approximation for 
the SNR. 


b. Obtained by using the 
Biot-Savert law and 
evaluating d8/^Z 


Z=0 


cont inued 

Legend for the table 4*2 j 
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E 

E 


1 


coil 


w 


1 


* 

2 

A Y 


®1XI 


9s 

6s 

R 

6 

0 

n 

r 

T 

Q 


■p 

^R 


£!tZ 


1. 


w 


E 


Sample characteristics 
®oil characteristics 
ojperatijag frequency ( oil^) 

•Signal decay time 
Sample volume 
(Soil teaiperature 
coil resistance 

aiesistance for sample induction loss 
asesistance for sample dielectric loss 
Sample radius 
example resistivity 
Sample skin depth 

field due to unit current in the coil 

-coil resistivity 

length of the coil 

coil radius 

s-kin depth in copper 

.angle substended by saddle at the axis of the coil 

number of turns in the coil 

radius of the wire 

dielectric loss factor 

coil quality factor 

coil distributed capacitance 

tuning capacitance 

pulse duration 

pulse repetition interval 

de ringing time 

object eleaient size 

observation window- 

current required in the gradient coils 
geometric factor for the gradient colis 
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CASE S!HJDY 

Tg = 310°K y/2n= 42.576 lEz /Tesla. 

T = 310°E; 

.■*1, = 3*25x10 ^ A m ^ (for protons in water) 
f Q = 1 .69x1 o'”® oiims. m 
= 0.85 ohms m (lOOmM) 

EaCl (Saline solution) 


Proton Imaging 


!YSTEM 

5 PECIFI CATION 

IMAGING SYSTEM EOR 

MJITI-PHASB 

BIO-SAMPISS 

CTLINICili 

IMAGING SYSTEM 


50-500m sec 

50 - 500 m sec 

^2 

Sample size 

50 psec- 500 psec 

1 cmx3cmx1 cm 

500psec-50msec [ I3] 

20cm2:20cmx25cm 

Dotal sample 
7olume Yg 

3x10"^ m^ 

O.Olm^ 

Eypical object 

contrast 

< 2 > 9 (r)>/< 9 (r)> 

2 - 5 f. 

2 - 5 ?^ 

,T2Variat~ 

ions 

ICP^ 

lof. 

Resolution 

required 

-v 5mm 

^ 5 mat 

Magnet 

available 

Iron Magnet (gap); 

5 cmx( Ti{ 2^) ^cm^ ) 

Saturation current 
^ 50 A 

Y^ = 0.02 

Eor a saddle shaped coH 

Y = 2 Y , gives uniform field 
wS-thifi 20 < over the sample [29] 

EP Probe 

Soleno id;Ri= 2 cm 
li= 4 cm 

Saddle ;R= 2 cm 
shape L= 4 cm 

0=120° 

saddle shaped ; R=l 5 cm 

Ii= 30 cm 

0=1 20° 
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amjnM umisa stsxbm foh 

SFFCIFICAIION M[JlTI-mSB 

- .. 3JstSimm 

'^INIOAIi 

IHACINa S7STBM 

Llue?r Oradieut 
aol3io‘vable 
( IDjrpieal) 

©•05- 
1 *00 /on 

within 109 ^ linear lt 7 

0«1 

-0 •20/om 

within 109 ( linearity 

Choice of field 


0«699HHz 

So 

f o( *50MHz( Saddle) 

1 

at f^a15MHz 


at f^ « 15MHz 

0 


-1 .5 

*0 sol 

17.21 


(1®) «'.0J6 

^0 saddle 

»8,Xf i 50HH!i 



20 MWt 


6g * 12 on 

1 9om 


I(15Maz,0.35 a!eala)/M 35 A 

- 


s 0*35 Tesla 

H s0.t4 Tesla 


In medical tm^l^ [9]> opti- 
mum field Is 3"4EB* due to 
lower 9 . 


Sim 

Sol. 

Saddle 




1 ,096x10“^ 

1 •151x10’^ 

3 . 24 x 10 "^° 


3?/% 

4.39acl0"® 

9-9x10“^ 

1 .9x10"^^ 



73-2dB 

67.7aB 

19dB 



(for T 2*2 .27 xl o"*^seo) 

(for T2ss500viseo) 



58.6dB 

52.22dB 

80dB 

■ ■ i 


^Vssl xl 0*^m^ 

(for T2*50ii8eo) 

(for Tg a 6.3vis®o.'^*=0-0l®^) 

Signal leweal 

7.150niiT 

11 .SnjiT 

0 .8411 wY 


(at taO) 

at a. 7 a 3rl O 


• 
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continued 




SI^GlPICAIIQlir 

Sesolut Ion 
a ^ 

I 

II!^= 500 iiseo 

X^sSOmseo 


<n 2>’/2 




SYSTEM FOB 


NDLTX-mSS 

..-jmTfiMSJiSS 


CLINICAL IMACINQ 
SYSTEM 


2.6xlO"^°m? 5.5xlO“^°m5 7.82x10"^n? 

2.6xl0“^^ii? 5.5xlO*“'^n? 7.82x1 o“^b? 

m ■ ■ 

• ^/n/( 1/T„) la a first order approxiaation, 
O' O' ' w' 2 .g / tp 

and typioally ^ 

0.01 In^ ohms 0.117n^ohins 1.26n^ ohms 
6.94npV/yBz 22.4np7/Vfiz 74p7/Vfia 


Pulsed Power. - 
Requ Irenents ( V^p) 

Broadband 


<100V pulsed rma 
(TpSslOiifleo, 0 ssl 8 O®) 


Without Gradient 25 OW 


(T *5jiseo, 9f*180°) 

Jr 


6 . 5 -low 


< 100 W 

(T^ss 500 (iseo, i^» 180 ®) 

» Jr 

*^ 6 mr 

■ • 

(Tp* 50 visec, 0 » 9 O®) 

.v 25 W 

• * 
.'j(Tpw 500 paeo, 0 * 90 ®) 

^IKW 

■ ■ (Tp«50|»seo., 0«9O®) 

/vl . 45 - 6 KW 

(Tpsl 00 - 50 viaeo, 0 » 18 O®) 

fM 350 W- 1 . 45 KW 

(T_*al00-50|iseo, 0=90®) 

Jbr 

0 . 4-1 OW 

for Tp/T;^ /X/ 10*^ 


for Tj^a/ 200^1890 



cont iaued 
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SYSTEM BTAGHTG SYSTEM POR 

SPSCIPICATIOM MJLTI-HIASE 

BIO-SAMPLBS 


Signal aquisition 
time 

A* 

Saaipling rate P 


Sampling 
interval T 


s 


Aquisition tire 
^a 


^700-3500 Hz 

-1 

1 


^ 2x10 
> 2x1 Ox 


27iT^ 


=70 KHz 
(T2 = 50psec) 

A^1 5psec 


6.57x10”^ sec 


2 1 / 2 ' 

(without averaging) (for 7^/<n > ' 


Maximum number 
of samples 


tv 40 


CLIHIOAI IMA.GIHG 
SYSTEM 

12.7 - 25 KHz 

1 .25 - 0.63x10”^ sec 
> 2x^^'* = 400 KHz 

2 .5psec 
lOOpsec 

^40 



HARDWARE AND SYSTEM SPECIFICATION 
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CfHAPTBR 5 


1 PULSED rf TRMSCBIVBR 

5.1 'TRANSKSCBIVER RSQUISEMERTS ; 

!I?h.e specifications for the design of the rf -coupling 

l. Out® corisicJeredL hare 

the transmitter and the receive^ She MR jmaging applications 
are discussed in chapter 4* Pahles4*5 and4»4 present a case 
study of an MR imager for multiphase hio -specimen, and the 
parameters to he specified on the transceiver* !I?he present 
chapter discusses the design of an rf coupling circuit and trens— 
ceiver Conf igurat ions for the application. 

The basic constituents on the rf transmitter-receiver 
are listed below : 

Transmitter 

1. Impedance matching with load (the sample coil 
resonant circuit), 

2 . High pulsed power from the transmitter amplifier, 

3. Small post-transm.it pulse deringing tine 

4 . Short power amplifier rise tine (or large band-width), 

5. Modulating signal design for selective pulse excitations 

6. Transmit switch offering low ’on ’attenuation and high 
»off * isolation, to decouple transmitter noise 
from the sample coil during reception of PID. 
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Receiver 

1* Low Rolse Figure for the preamplifier, 

2* Efficient noise, and signal power matching of the 

source (the sample coil) with the receiver preamplifier, 
5* Quadrature demodulator for the recovery of the base 
band signal (with anti-aliasing filter), 

4* xf switches with high isolation to protect the 
receiver from high power rf pulses and low 'on» 
attenuation, and shot noise contribution, during 
signal reception. 

The similarity with the Radar Context; may be noted. 
Figur^.4 illustrates a typical excitation - reception 
seq.uence in imaging, with the various timings involved therein. 
The i!able5»1»lists the specifications for the system under 
development. The deringing time is estimated for the worst 
case values of T 2 ^ and T^* The resistance (R), inductance (1), 
quality factor (Q^ max*.) and the distributed capacitance 0^ of 
the coil are given along wllhpsuk and average current (I and 

XT 

^ayg^ and the rf field (B^) requirements for T^ 5p.Sec and 

0 V2* 

The specifications for the BNI-series broadband power 
amplif.iers and the Matec series rf preamplifiers are assumed 
in the design of the coupling circuit given in the following 
section. 
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5.2 rf COUPLING; CIRCUIT 

The design is based on the following observations ; 

1* The picknp coil (same as the transmitter coil, for 
s in^e coil probes) is an inductive source (or load) with 
an inductance given by ls=Cn , where C is a constant detemined 
by the geometry and the dimension of the coil, and n is the 
number of turns in the coil. Here n is a variable, within 
physically realizable limits^ 

2. The received signal and the excitations are narrow 
band, so that tuned circuits may be eci ployed* 

5* The excitation rf field is propoitional to the 
current in the coil, which is maximum in series resonance 
Circuits for a given external driving voltage. 

Consequently, the circuit design for series resonant 
circuit is considered.in this thesis. A starting point is the 
circuit shown in Pigure5‘t, due to Clarh and McNeill [30]*. 
Three sets of independently tuned circuits are employed as 
shown* Cj and Cjj are chosen to meet the transmitter 
requirements while is obtained from noise matching 

Considerations. The switches and Sg are the transmitter- 
receiver switches, closed during excitation and open during 
reception. The transmitter coupling design is obtained by 
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maxiniziag the field with the sample coil I 2 , within the 
constraint of An assumption, that only half the energy 

in the coil is stored within its confines, is made to relate 
with {Figure 5*^) • This is not valid for saddle shaped 
CO ilSj where the field concentration is near the conductors* 

The results of the calculations jare indicated in Table 5*2 . 
The choice of I 12 is made based on aialLu n number of turn® roiuirod 
on a solenoid for nominal field homogeneity within the coil* 

The field ache ived is 0.75 B* whore B. is the 

maximum field for a given current 1,^1 *• The design indicates 
large power dissipation (**> 3k¥ pulsed mns) in the external 
resistor R-g » 21 ohms, employed to reduce Q^. ilfhe impedance, 
mat’dh, = ■ 58*5 “ 327*9 ohms, is also poor. As such this 
transmitter coupling circuit configuration is not well suited 
for imaging where power reiuiroments are large. 

The receiver' coupling configuration R^, is optimum 

for noise matching of narrow band resistive signal sources with 
the preamplifier (T is hay Fetzer The design is Carried 

out assum.ing an e^ (equivalent input series noise voltage 
source) and ^ (the equivalent input shunt noise current source). 
The SFR, optimum ratio (1^/C^ for minimum SIR) and the matched 
noise f igure (H^ are expressed, as ; 
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sm = 


4WEg + Eg\2 (1 + 1,/Hg‘^C) + Eg^OA 




(lg/C^)Qp^ - Rg /®I1 /^ - ^ ■» (^^*^■5) - 

^^^matched 

Por typical PST input stages, 10-25nV/Y-^^ and i^ /■..» 1 pA/ y^^ » 
This yields a noise figure of '^6dB, (L^/C^) ^ 25500, while . 

Ii^ and 0^ are 1 *^pH and 60-^0 pf respect ively. As such given 
preamplifier, a Imowledge of and i^ are required to 
obtain noise matching. The Matec 615 tuned receiver available 
guaranties 6dB noise figure at 50 ohms input impedance. The 
variation of the UP around the optimum ratio is small for low 
noise amplifiers (Uetzer [31]),. As a result small mismatches 
do not affect the noise perfoimance of the amplifier. 

An alternative coupling circuits is shown in figure 5-2, 

The impedance matching is obtained by the capacitor while 

Q control is effected through the use of a resistor 10 ohms. 

o 
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The residual reactance 3^ = tuned out with the 


inductance through impedance matching. The series circuit 

•J 

£ 2 # I' 2 , C 2 is detuned by a factor £> = ■WI 2 ~ , to obtain 

2 .4 

matching with 0^* To keep the residual reactance at lO ohms 

n 


A must be equal to 5« Consequently for receiver tuning the noise 


matching network 1^^. is chosen around the optimum to tune out 
the source reactance jA . This circuit shows proper impedance 


matching and lower power requirements, though it has the dis 
advantage of a complicated tuning procedure (figure 5*2). 


Pigure5«5#hidicat9S a circuit wherein the deringing, after 
the cessation of the pulse, is reduced using a pulse quench ■ 
circuit with the switch Sq* The conf iguration during reception 
is indicated in f igure 5 The sources matched to V4 cable 
by the split capacitance pair C 2 , (used for transmitter 
matching) . The cable with the crossed diode pair D is chosen 
to protect the receiver from the high power rf pulses* The 
transformer of turns ratio 1 :n is used to effect noise matching 
with the preamplifier. 


The design values are given in Table 5 •2* While the 
circuit parameters I 2 , O 2 , and u "the turns ratio of the 
transformer are realizahle, and the pulsed power requirement 
¥p < 500¥, the reactive m.ismatch in is 10 ohm.s, which is 
a significant value. This is constraingdby the problem of 
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pulse <laench.ing* The reiuirement here is that the *on’ resis- 
tance of the luench switch tnust he STOallei' than the capacitive 
impedance* Ideally is kept large enough to obtain a compro- 
mise between pulse de ringing and reactance mism.atch* 


This problem is alleviated In a similar circuit by- 
Kisman and Jknstrong [ 32] | "by use of TIN diodes as a part of 
the resonant circuit to effect quenching. However the rf 
power dissipation and shot noise considerations ( "^GOpA/^^^) 

make its scope limited in the content. Hoult [ 53], presents 
a scheme, with TIN diodes turned off during reception* This 
Circuit involves the use of four quarter wavelength sections 
( X /4 rj 5m at 15 MHz), which is undesirable at the frequency 
of interest* 


To sum up, the choice of coupling circuits with series 
resonance circuits involves a trade off between power matching 
pulse deringing and receiver noise figure requirements. . 'Where 
the eristance of good switches are assumed the scheme in 
figure 5*5 1 .is feasible, else the first scheme may be resorted 
to* 


A wide band ImpGcbnce transformer, reducing the source 
impedance of the amplifier (E.^'^50 ohms) to a lower value, 
would introduce flexibility in the choice of the matching 
capacitance C^* The matching need now be done for the 
receiver circuit only. 
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TiSLB 5.1 


SYSTEM SPECIE ICAT lOE 


PARAMBTEE 


^ 2 '' (for 0.05 G/cm) 
T^ 


T 


1 


R 

1 

Q 


c max 

f 

'd 

'■avg 


.SQEYMOID 

0.01 In^XL 
0.027n^tPS 
200 
^ lOpf 

105.8/n A peak 

(,W) l06mA/n. 

avg 

^“5 


VAIUE 

2.27 z 10 ^ sec 

50jaSec ~ 5niSec 
5mSec - 500mBec. 

A) ImSec 

"6 

— ^ 0.5pSec- 

s1ri)le~^7ap¥d' 


. 5.QzlQ 

^ 100 


0.1 76n ja. 

0.099n^pH 

32 

< iopf 

70/n A peak 
fY2-)70/n . avg. 
-5. 


2.22x10 n tesla 3*67x10 n Tesla 


(B^ is calculated for 0 = 90*^, T = 5 sec) . 




cm cotr A a Tr<A»iSM)TT£ft oK; s, on’ ^ % on 

Ti^ANiSMiTTcR Qfv^ *> 0('f- 

— C*- 2 .^:-) “ C ^ "3 ) ' ^'C 

ifi THt SA>APUE. <ra(L. 

Ff^ORE ^ i A SE F<»\£S fsC <IooPut\K^ 

ORCOIT CctAR< AND V^-cMeiU. f.'i'i')) 




t 


ciRcoiT 8 •; a^l-^ - J_/j_ , l.^ ,v o 

'-,rr,> C2J 

- _L- =Q 

R: ^G' SPon.INfi RFSt-STANC.t: 
s 


l=}GOVRS 5'>1 RfcSoK<AJ^iT ClfNCOV"' VNnH voqsN&k 
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RE'CE.PTtOt-i 

% Fotse ^WtTCH 

2o = ^ j_ A =: f tot, -1_ \ 

wG-fts. 2*>Cw» c^C2J 

SPUT CAPAOTAf^Ct HATOHrCj 

FitaORS- '5*3 COUPi-lf'Jii CiRCVlT UtSJrsJiS 

SttTSOlO P/K t><O0E SJ^ITCH 

-**'f7p l-sr- 



FC >•>'•'> TU'^N'b' 


P>6oPC 5./.» Topical e xciTAnorNS SiauE^’JCES 
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liBLE 5.2 



CIRCUIT 

VALUE ROR A,B, 

and C. 

C ircu it 
Elements 

A . 

B 

0 

^2* ^2 ^^2 

0.4 IliH 1 12..5pf 

0.4 llxH, 120pf 0.4 IpH 120pf 

-l/V ^ c- -lii. 

A =5 

Rl,I^ ,C^Rg 

^ , 0 ,' 5 \ i 3 .. 579pf 
,21/1- 

- 

- 


- 

1 .ItiH, 10^2. 

I05pf 

0,0, Ikpf 


in-,r.7~2|iH, 

60-^0pf 

in., 1 .7-2pH, 
40-60pf 

%pt = 22.5 

^i, 

38.5-327.99^ 

50 Xn . 

50 “jio.ort- 

% 

-*6dB 

^ SdB 

Aj6dB 

¥ 

P 

/^3.1k¥ 

^ 1 .5k¥ 

^ 60--250V 

'Wg 

^ 3.1W 

^ 1 .5¥ 

- 


6 

AlO 

Q 


>. 50 

>50 

a; 

50 


f: 

= 50 t' - 2 7 l -< 15x10 rad/sec 
= 0.5iiSec, n -= 6 
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TABIfS 5*2 (continued) 


EQUATIONS FOR THE CIROUITS A,B and C. 


\ < WL 2 Qg, iL 1 0 R^ 

WT = 2 Q„ 


®- a s ''U'i+i' 2 ) 


I't = . (H1+K2’ = 




B, > 2/5 I,/C,={E2+II,) /l^ 

\ = ^ 


c_ = 


A 

'm T^^Sg + Rp“ 


A = WL2 ~ 1 /'”’°2 

I'3) and (©2, l2» V indepen- 

dently tuned at w. 


Z. = 


- 1 * 


WC^ R2 V/ 
A 

°m ~ 'KT R^ X“ 


"0 ^'2 


WC ^0 

m 


^ = ’*2 " sj: 


(c , C2, Iig) tuned. 



CHAPTER 6 


STOMART 


6.1 DISCCSSIOH 

lUhe present work discusses the MR amaging method 
from the point of view < of system design* The ohservations 
ensu-ins ' frcan a mathematical model of the imaging system* 
and the system design considerations are discussed. 

The model based on Bloch equations and rotating frame 
arguments show that, 

1 * iKie latter is a fairly useful tool in modelling 
different modalities in MR imaging. 

2 * The response (namely the magnetic moment m_ (t)) 
is. linearly related to the system identification function 
) L(r^), while the aperture function is non-linearly 
related to the rf excitations-. 

5. The nonlinearity is found in the Bloch equations 
in the co-efficient B^(t). The solution for a rectangular" 
excitation is possible in closed form, while for other pulses, it 
is obtained numerically. 

4* The choice of pulses for signal design (in selective 
excitations) can be made by approximating the system by the 
linear teim M ^^^(t) (section 3-4)* The contribution of 
the linear term to the overall response is s ignif icant and 
is 50*^ for rectangular pulse excitations (90^ flip). 


Ill 


3 » Under white gauss ian assumption (on the quadrature 
noise components), the estinate of the spin density function, 
for broad hand excitations is 

(w) = yCw) + ITCw) where U(w)isawhite gaussian process. 

The signal to noise ratio for a Zeugmat ©graphic experiment is 
related to (the signal energy-to -noise spectral density), 

as in finite time signals* The ratio is given by 




o 



for the case of a strong linear gradient and nominal contrast* 
This is similar to the SNR for a single pulse MR experiment. 


The system design considerations indicates that, 

1*. The hierarchy is in three levels, namely the 
application, object specification and syst^x parameter speci- 
fication respectively,. This schematic is useful in design as 
well as in evaluating existing systems for a given task. 

, 2. OJhe distinction between the detection and display 

stages in imaging, enables the evaluation of the basic image 
quality (in terns of resolution and contrast) determined by 
the system hardware . 
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Tiie resolution and contrast in the detection stage, 
are limited by noise. /While in reality the tifo parameters 
are interrelated^ a first order estimate is carried.* -out by 
considering point pair resolution and noise (in terms of SWR) 
in Isolation. . 

, 4<* IPw-o basic parameters associated -with the detection 
stage are 

a* Signal to Ho ise Ratio 

b. !Phe choice of operating frequency. 

fhe maximum SHR obtainable in an imaging experiment has the 
following proportionalities 

a. Equilibrium magnetization ( ct M^) 

2 ' ' ’ 

b. Volume of the sample ( a V_ ) 

s 

6. The signal decay time (a T^. ) 
d. The frequency w^ {« 

e* Signal averaging gain (a 

when the only non— ideality considered is the skin effect in the 

2 

coil, which affects the figure of merit /^c* 

At high frequencies where inductive losses in the 

sample are predominant , the frequency dependence becomes 

w The SIR advantage in frequency, and the hardware complexity 

o 
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(in teiws of the magnet) define the optimum field range in 
imaging. S^or typical systems it is 3“403-^ 

5* Ihe received signal level is of the order of a 
microvolt^ -t/hile . the m in imito transmitted power level required 
is '^1' (20/yT^)^ -£ 2 .^ pulsed rms * Ihis may be more than 

■ ik¥ in typical systems* Considering the dissipation in 
other parts of the system j the total power required from 
the amplifier can be < tOk¥« 

For the case of imaging systems for multiphase bio - 
specimen, 

1. Con'vent ional pulsed, wideline NMR spectrometers 
may be modified for imaging. 

, 2k A m.il3.imeter resolution is possible with a gradient 
of 0.5 G/cm, assuming a 200 Hz frequency resolution and 
signal averaging* 

3* A feasible modality is the projection reconstruction 
technique, due to the simplicity of data aquisition. 

4* The coupling circuits, with single coil probe, 
show stringent requirements in terms of impedance matching, 
power dissipation and post pulse deringing tlm.e* For large 
power requirement, adequate receiver isolation and quenching 
is called for* 
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5* Ihe use of PUT diodes, for series transceiver gate 
and q.uen.ch switches, is not desirable due to its large shot 
noise contribution lOOpA/yjj^) • Shunt switches are cumber- 
some for lower frequencies ( I 5 MH 2 ), due to a number of 

^4 sections required* As such, where space is not a premium, 
crossed coils offer advantage in teims of convenient transceiver 
in pi em ent at ion » 

6.2 SUGCtESPIOITS FOE future WORK 

1# Evolution of the imaging equation covering situations 
such as the chemical shifts, which is especially- relevant in 
study of in -vivo metabolism, 

2* Identification of selective pulses using the linear 
assumption in Bloch equations, and numerical solutions. 

3- Study of second order effects, such as gradient 
non-linearity, field fluctuations, and rf non^^nifo unity over 
the sample, in system identification. 

4 . Derivation of a contrast-resolut ion-SllR relationship 
for different imaging tasks, to -''eteimine realistic bounds 

on the performance of the system. 

■ % 

5 . Designs of coi3.vS, for the magnetic fields (Hq!< 6E&) 
and B^(t) , from homogeneity and SNR ( a ^/R^) , and linearity 
(in field gradient) cons idorat ions . 


6- Design of a sin^e coil probe, for imaging, 
catering to the impedance matching and deringing requirements. 

7*. A realistic analysis of the optimum field ranges, 
by simulation of different physiological samples. 


EBFSEBITGBS 


1# lauterbur* P, 'Image Formation by Induced Local 

Interactions r examples employing UMR', Hature, Yol.242, 
No. 5394, PP 190-191, Ifercb 1973* 

2, Damadian* R, 'Tumor Detection by MR', Science, Y0I.I7I, 

Ho .3976, PP 1151-1153, March 1971. 

3. Moss. A, 'Imaging Modalities', lEBB Trans., Medical 
Imaging, MI-1, Ho. 4, PP 215-219, December 1982. 

3. Macovski. A, 'Selective Rejection Imaging : Applications 
to Radiography and HMR', IBBl Trans., Ifedical Imaging, 
MI-1, No.l, PP 42-47, July 1982. 

5. Ter-Poggosian. M.M,, 'Bngineering Research in Medical 
Imaging : A Synergistic Relationship', IEEE Trans., 

Medical Imaging, MI-1, Ho. 4, PP 205-209, December 1982. 

6. Cho Z.H. et.al., 'FT-HIIR Tomographic Imaging', Proc. 

IBBB, V0I.70, Ho. 10, PP 1152-1173, October 1982. 

7. Kaufman et.al., 'Physical Basis of HM Imaging’, IBBB 
Trans., Huclear Science, HS-29, Ho .3, p 1252, Juno 1982. 

8. Scudder H.J , 'Introduction to Computer Aided Tomography’, 
Proc. rSBB, Y0I.66, Ho .6 , pp 628--637, June 1978. 

9. Kaufman. L and Crooks L.B , 'Realistic Expectations for 
the Hear Term Development of Clinical HMR Imaging', IBBB 
Trans., Medical Imaging, MI -2, Ho .2, pp 57-65, June 1983* 



117 


lO* Andrew E.R , ’Euclear I^Iagnetic Resonance*, Cambriclge 
Monographs on Physics, Univ. Press 1956. 

11 . Ernst R.R. and Anderson ¥.A , ’Application of Pourier 
Transform Spectroscopy to Magnetic Resonance’, Review 
of Sclaatlfid Inffbruinsnts^RSI--^?, Ro.l, pp 93~102, 

January 1 966 • 

12, Hoult D,I , ’ Zeugmatography : Criticism of the Concept 

of Selective Pulses in the Presence of a Field Gradient’, 
Journal of Magnetic Resonance , JMR-26, pp 165 “167, 1977« 

13 • locher P*R , ’Proton RMR Tomography’, Philips Tech. 

Review, Vol. 4i, Ho.5, pp 73~88, 1983“84« 

14* .Mansfield P. et.al,, ’Fast Scan Proton Density Imaging 
hy RMR*, Journal of Physics, 3. Sci “ Instrum Yol.9, 
pp 271 “278, April 1976. , 

15- Ifeudsley A. A , ’Multiple Lino Scanning, Spin Density 
Imaging’, Journal of Magnetic Rosonance, JMR“4l , Uo-I, 
pp 1 12-126, October 1980, 

16. Kumar A. et.al., ’RMR Fourier Zeugmatography’, Journal 

of Magnetic Resonance, JMR-18, Ko*1, pp 69“83, April 1875* 

17* Hoult. D.I , ’The Solution of the Bloch Equation in 
the Presence of a Varying B-) Field - An Approach of 
Selective Pulse Analysis’, Journal of Magnetic Resonance, 
.JMR-35, Ho.1, pp 69“86, July 1979- 


118 


18. ffetnsfield et*al., ’Solc-ctivo Pulses in UMR Imaging 

A Eoply to Criticism’, Joumal of Magnetic Rosonance, 
3rMR-33, No *2^ pp 261-274, February 1979* 

19* Oaprihan,. A, ’Effect of Amplitude Modula,tion on Seloctivo 
Excitations In NMR Imaging’, IEEE Irans, Medical Imaging, 
MI-2, No. 4, pp 169-175, December 1983. 

20* Hinshaw ¥.S* and Lent A.H., ’Bloch Equation and tho 

Imaging Equation*, IBEE Proc., Yol-71, No .3, pp 338-350, 
March 1983^ 

21* Silver et.al., ’Selective Pulse Creation by Inverse 
Solution of the Bloch - Riccati - Equation’, Proc. 

2 nd Meeting, Society of Magnetic Rosonanco in Medicine, 
1983. 

22. Hoult D.I. and Richards R.E., ’The SNR of an NMR 
Experiment’, Journal of Magnetic Ecsonanco, JMR- 24 , 

N 0 .I, pp 71-85, Octobor 1976. 

23. Hoult D.I. et.al., ’Bloctro magnet For Nuclear Magnetic 
Resonance Imaging’, Review of Scientific Instruments, 
RSI~52, No. 9, pp 1342 - 1351 , Septenbor 1981. 

24. Hagnor R^F. and Brown D.G., ’Overview of a Unified 

SNR Analysis of Modical Imaging Systems’. IEEE Trans., 
Medical Imaging, MI-1, N0.I, pp 210-214, December 1982 * 


119 


25. House ¥«7», ’HMR Microscropy’ , IEEE Trans., Huclear 

ScioncG, HS"-31 1 *1 » PP 570-576 , February 1 984* 

26. Hoult E.I* and Lautorbur ?«G., 'IHig Sensitivity of tbc 
Zougnatograpbic Exporinont Involving Human Samples ’ , 

Journal of Magnetic Sosonanco, JMR-34, Ho .2, pp 425-453, 

May 1979. 

27 . Honney E. cd., ’Radio Enginooring Handbook’, McG-rav Hill, 
5th odi-feiotij 1959* 

28. G-inzborg et.al., ’Optimum Goemotry of Saddle Shaped 
Coils for Gonorating Uniform Magnetic Fields’, Eeviev 
of Scientific Instruments, RSI*“4l , H 0 .I, pp 122-123, 

January 1970. 

29. Bernardo Jr. M.L. et.al., ’Radiofrequency Coil Designs- 
for Nuclear Magnetic Resonance Zeugmatographic Imaging*, 
Proc. International Workshop on Physics and Engg. in 
Medical Imaging, March 1982. 

30 . Clark G.¥. and McNeill J.A., ’Single Coil Series Resonant 
Circuit for Pulsed Nuclear Resonance’, Review of Scientific 
Instruments, RSI-44, No .7, pp 844"857, July 1973* 

Netzer Y., ’The Design of low Noise Amplifiers’, Proc. 
lEBB, Vol. 69, N 0 . 6 , pp 728-741, June 1981. 


31 


120 


52 • Kisman E.B. and Armstrong E.L., 'Coupling Scheme and 
Irobe Damper for Dulsed Nuclear Magnetic Resonance 
Single Coil Probe’, Review of Scientific Instruments, 
R3I-45, No. 9, pp 1159’*1162, September 1974* 

33* Hoult D.I. and Richards R.I1», 'An Ultra High Prequency 
Receiver Scheme’, Journal of Magnetic Resonance, 

JMR-22, No.1, pp 561-^73, April 1976, 

34* Brunner, P» and Ernst R.R, 'Sensitivity and Performance 
Time in NMR Imaging’, Journal of Magnetic Resonance, 
JMR-33, No.1, pp 83~106, January 1979* 


APPENDIX A 


ROTATING FRAME ARGUMENT (Rectangular pulse ezcitation) 

Some details of the equations (‘^•7 "'‘^•9), of this 
thesis are provided here. The coordinates are labelled 
X',T' and Z’* The tvo iateracting magnetic fields in the 
rotating fraiae are ; 

1 . ah, the residual field in the rotating frame 
(=G X, for a gradient along the X-axis ) . 

2. Bj(t), the rf modulating -waveform (=B-] , a constant 
for 0 < t < T for a rectangular pulse.) . 


y 

Figure A.1 Rotating Frame at Point x 



For figure A.1 , 


L 

lores' = ((V)^ 

G = tan”^ iB^( t)VG-j.(3c)) , the tilt angle 

^r = l^Jesi ^p » flip angle. 
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For a general initial condition, i.e., 0^ = 0|^_^(0o), 
the projections of the magnetization along the X^,T’ and Z’ 
axes are, 

= M^(x) cos ^^in© - sin cos (0^+0^ ) cos©"} 

= Mq(x) I sin§ sin (0^+0^)^ 

a nd ( x) s in ^ cos ( 0^ 0^ ) s in© 

cos'^ cos© ^ 

where ^ is the semiconical angle (or the angle hetween the 
vectors and M^(r,t), and 0^ is the flip angle). 

The projection or the aperture function F(x), expresse-d in 

terais of M , and M , is 

x y 


pTj:) = -4-.- (IL.-jM -,) 

M^(x) ^ y 

=: Sin © cos © (1 “ COS0) - j sin © sin 0 
for the condition 0^ = 0 and ^ =©» 

Given >> G^x(the broad band excitation), and 0^ = ti/B 



12^ 


and 


P^' = 0 


= \(^) at t=T^ 


(A.l) 


while for a5 


P^,(x) = sin e cos©’ (1-cos{^ 

‘ 3 . \aA 


and 


F ', (x) = sin ©• sin 0 

V 


= B, V Tj, sine (YfB2+(e‘z)'^ TJ 


P' 
(A.2) 


The Bloch equations and the geometrical arguments nay 
he related hy the fact that, the former reflect the dynamics 
in the X*,Y’ and Z* frame, while the latter treats the p 2 x>blero 
in a frame tilted with respect to the X’,Y’ and axes. 

The -vector B^„„ is along one of the axes of this frame. 

-L ©S 

The trajectory of the magnetization can he treated as a free 
precession about the direction . Further the path may 
he referred to the rotating frame X',Y’, and Z’, through a 
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transformation of coorrdinates • Ihus the solution of the 
Bloch equations is equivalent to a that of a solid geometry 
problem* 

For example, for^JiH « B. , M , = -j M (x) sin 0, as 

' max j u - 

may be inferred from the geometric arguments. The same result 
is obtained by solving the Bloch equations (for ji::>H=0) . The 
perturbation term M^^^(ln equation takes tho form 

n ^ 

^ M^(x) for n, odd and n^1 - 

u# 

This is the Taylor expemsion for My. 7 . ¥hen nH^O, the per- 
turbation series -would converge to equation (A. 2)* 



APPUFDIX B 


IM&ING- BQUATIOBS FOR SOJE MODALITIES 


1 . EBCONSTEDCTIVE TSCHRIQUES 


a. 3~D Plane integral projections 


in(t) 


= / 


e ^ ® dV 




f(r") 


s 


"b* 2~D lane integral projections 


{B.l} 


m(t) = K jiZ J ?( r ) e dA 


A 


s 


('B.2) 


(tlie actuaA reconstruction is visualized in polar coordi- 
nates) . 

2. ROI-RSCONSTRUCTITE TECHNIQUES 
a. Sensitive point scanning 

/ - jY( A(x)s-4-B(y)4-c(z)z) cos vrt 

f ( 50 e 47 

(B-") 

(The frequency modulation may he noted) 
h- Line scanning 


7/ 


dY^^xt 


m(t) = JJJ I' (x.,y,z)F^(z,z^)F^y(y,y^)e ^ dzdydz 

'(B-4) 


X y 2 
s^s s 


(Here the readout gradient is in the x-direction) 



E® Fourier Zeugmatography : 


ai(t) 


f 3y(G rb -!-G^ +(?2zt ) 

= J f(x,y,z) e 3c X y Z dxdydz 


V, 


s 


{B.5) 


where are 0 

RSCOVIRT SEQUENCES 
a. Saturation recovery ; 


m(t) 



(z) 

(1-e ^ ^ 


_3^/2 

)e ^ e 


dx 


whe re , 

i) <'E^(z)> 

^zy 

and iii) Spin lattice relazation has a exponent oal 
behaviour at point; z« (B.f) 


b . Inversion Recovery ; 

At tiia» t = Tj after the 90*^ pulse, a broadband 
180® pulse is applied, 



1 27 


and at t = T, 




-^/T^(z) 


-T/TiCx) 

+ M^(z) (1-« V ) 


. / \ r — T /T { x) f 

or M^(T,x) =M^(z) [ e -^/^r^^C-a+e ’ 


) + 1 ] (B. 7 ) 



E£EATA 


1 . 'Plank* appearing in several places in this work should 
read 'Planck*. 

2. 'Iiannour* appearing in several places in this work, should 
read 'Iiarmor*. 

3* ’Atuisit ion * appearing la several places in this work 
should read 'Acquisition*^ 

4* P‘77i -line 12: 'T.* should read 'T ' . 

5. p.79, line 20 J | 2 

1 *w ' should read 'w * « 
and p.tl2, line 21 1 

6. p.94, second row ; under flinical Imaging System', 

*f^ = 15 MHz* should read *f^=6MHz*. 

7* p.110, line 22 : '309&' should read '765&». 


